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Basics  of  Radar  Polarimetry 

Abstract  A  comprehensive  overview  of  the  basic  principles  of  radar  polarimetry  is  presented.  The  relevant 
fundamental  field  equations  are  first  provided.  The  importance  of  the  propagation  and  scattering  behavior  in 
various  frequency  bands,  the  electrodynamic  foundations  such  as  Maxwell’s  equations,  the  Helmholtz  vector 
wave  equation  and  especially  the  fundamental  laws  of  polarization  will  first  be  introduced:  The  fundamental 
terms  which  represent  the  polarization  state  will  be  introduced,  defined  and  explained.  Main  points  of  view 
are  the  polarization  Ellipse,  the  polarization  ratio,  the  Stokes  Parameter  and  the  Stokes  and  Jones  vector 
formalisms  as  well  as  its  presentation  on  the  Poincare  sphere  and  on  relevant  map  projections.  The  Polarization 
Fork  descriptor  and  the  associated  van  Zyl  polarimetric  power  density  and  Agrawal  polarimetric  phase 
correlation  signatures  will  be  introduced  also  in  order  to  make  understandable  the  polarization  state 
formulations  of  electromagnetic  waves  in  the  frequency  domain.  The  polarization  state  of  electromagnetic 
waves  under  scattering  conditions  i.e.  in  the  radar  case  will  be  described  by  matrix  formalisms.  Each 
scatterer  is  a  polarization  transformer;  under  normal  conditions  the  transformation  from  the  transmitted  wave 
vector  to  the  received  wave  vector  is  linear  and  this  behavior,  principally,  will  be  described  by  a  matrix 
called  scattering  matrix.  This  matrix  contains  all  the  information  about  the  scattering  process  and  the 
scatterer  itself.  The  different  relevant  matrices,  the  respective  terms  like  Jones  Matrix,  S-matrix,  Muller  M- 
matrix,  Kennaugh  K-matrix,  etc.  and  its  interconnections  will  be  defined  and  described  together  with  change 
of  polarization  bases  transformation  operators,  where  upon  the  optimal  (Characteristic)  polarization  states  are 
determined  for  the  coherent  and  partially  coherent  cases,  respectively.  The  lecture  is  concluded  with  a  set  of 
simple  examples. 

1.  Introduction:  A  Review  of  Polarimetry 

Radar  Polarimetry  {Polar,  polarization,  Metry:  measure)  is  the  science  of  acquiring,  processing  and 
analyzing  the  polarization  state  of  an  electromagnetic  field.  Radar  polarimetry  is  concerned  with  the 
utilization  of  polarimetry  in  radar  applications  as  reviewed  most  recently  in  Boerner  [1]  where  a  host  of 
pertinent  references  are  provided.  Although  polarimetry  has  a  long  history  which  reaches  back  to  the  1 8th 
century,  the  earliest  work  that  is  related  to  radar  dates  back  to  the  1940s.  In  1945  G.W.  Sinclair  introduced 
the  concept  of  the  scattering  matrix  as  a  descriptor  of  the  radar  cross  section  of  a  coherent  scatterer  [2],  [3]. 
In  the  late  1940s  and  the  early  1950s  major  pioneering  work  was  carried  out  by  E.M.  Kennaugh  [4,  5].  He 
formulated  a  backscatter  theory  based  on  the  eigenpolarizations  of  the  scattering  matrix  introducing  the 
concept  of  optimal  polarizations  by  implementing  the  concurrent  work  of  G.A.  Deschamps,  H.  Mueller,  and 
C.  Jones.  Work  continued  after  Kennaugh,  but  only  a  few  notable  contributions,  as  those  of  G.A.  Deschamps 
1951  [6],  C.D.  Graves  1956  [7],  and  J.R.  Copeland  1960  [8],  were  made  until  Huynen’s  studies  in  1970s. 
The  beginning  of  a  new  age  was  the  treatment  presented  by  J.R.  Huynen  in  his  doctoral  thesis  of  1970  [9], 
where  he  exploited  Kennaugh’s  optimal  polarization  concept  [5]  and  formulated  his  approach  to  target  radar 
phenomenology.  With  this  thesis,  a  renewed  interest  for  radar  polarimetry  was  raised.  However,  the  full 
potential  of  radar  polarimetry  was  never  fully  realized  until  the  early  1980s,  due  in  no  small  parts  to  the 
advanced  radar  device  technology  [10,  11].  Technological  problems  led  to  a  series  of  negative  conclusions 
in  the  1960s  and  1970s  about  the  practical  use  of  radar  systems  with  polarimetric  capability  [12].  Among  the 
major  contributions  of  the  1970s  and  1980s  are  those  of  W-M  Boerner  [13,  14,  15]  who  pointed  out  the 
importance  of  polarization  first  in  addressing  vector  electromagnetic  inverse  scattering  [13].  He  initiated  a 
critical  analysis  of  Kennaugh’s  and  Huynen’s  work  and  extended  Kennaugh’s  optimal  polarization  theory 
[16].  He  has  been  influential  in  causing  the  radar  community  to  recognize  the  need  of  polarimetry  in  remote 
sensing  applications.  A  detailed  overview  on  the  history  of  polarimetry  can  be  found  in  [13,  14,  15],  while  a 
historical  review  of  polarimetric  radar  technology  is  also  given  in  [13,  17,  18]. 
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Polarimetry  deals  with  the  full  vector  nature  of  polarized  (vector)  electromagnetic  waves  throughout  the 
frequency  spectrum  from  Ultra-Low-Frequencies  (ULF)  to  above  the  Far-Ultra-Violet  (FUV)  [19,  20]. 
Whenever  there  are  abrupt  or  gradual  changes  in  the  index  of  refraction  (or  permittivity,  magnetic 
permeability,  and  conductivity),  the  polarization  state  of  a  narrow  band  (single-frequency)  wave  is 
transformed,  and  the  electromagnetic  “vector  wave”  is  re-polarized.  When  the  wave  passes  through  a 
medium  of  changing  index  of  refraction,  or  when  it  strikes  an  object  such  as  a  radar  target  and/or  a  scattering 
surface  and  it  is  reflected;  then,  characteristic  information  about  the  reflectivity,  shape  and  orientation  of  the 
reflecting  body  can  be  obtained  by  implementing  ‘ polarization  control’  [10,  11].  The  complex  direction  of 
the  electric  field  vector,  in  general  describing  an  ellipse,  in  a  plane  transverse  to  propagation,  plays  an 
essential  role  in  the  interaction  of  electromagnetic  ‘ vector  waves  ’  with  material  bodies,  and  the  propagation 
medium  [21,  22,  13,  14,  16].  Whereas,  this  polarization  transformation  behavior,  expressed  in  terms  of  the 
“polarization  ellipse”  is  named  “ Ellipsometry  ”  in  Optical  Sensing  and  Imaging  [21,  23],  it  is  denoted  as 
“ Polarimetry ”  in  Radar,  Lidar/Ladar  and  SAR  Sensing  and  Imaging  [12,  14,  15,  19]  -  using  the  ancient 
Greek  meaning  of  “ measuring  orientation  and  object  shape”.  Thus,  ellipsometry  and  polarimetry  are 
concerned  with  the  control  of  the  coherent  polarization  properties  of  the  optical  and  radio  waves, 
respectively  [21,  19].  With  the  advent  of  optical  and  radar  polarization  phase  control  devices,  ellipsometry 
advanced  rapidly  during  the  Forties  (Mueller  and  Land  [24,  21])  with  the  associated  development  of 
mathematical  ellipsometry,  i.e.,  the  introduction  of  ‘ the  2x2  coherent  Jones  forward  scattering 
(propagation)  and  the  associated  4x4  average  power  density  Mueller  (Stokes)  propagation  matrices’  [21]; 
and  polarimetry  developed  independently  in  the  late  Forties  with  the  introduction  of  dual  polarized  antenna 
technology  (Sinclair,  Kennaugh,  et  al.  [2,  3,  4,  5]),  and  the  subsequent  formulation  of  ‘the  2x2  coherent 
Sinclair  radar  back-scattering  matrix  and  the  associated  4x4  Kennaugh  radar  back-scattering  power- 
density  matrix’,  as  summarized  in  detail  in  Boemer  et  al.  [19,  25].  Since  then,  ellipsometry’  and  polarimetry 
have  enjoyed  steep  advances;  and,  a  mathematically  coherent  polarization  matrix  formalism  is  in  the  process 
of  being  introduced  for  which  the  lexicographic  covariance  matrix  presentations  [26,  27]  of  signal  estimation 
theory  play  an  equally  important  role  in  ellipsometn’  as  well  as  polarimetry  [19].  Based  on  Kennaugh’s 
original  pioneering  work  on  discovering  the  properties  of  the  “Spinorial  Polarization  Fork”  concept  [4,  5], 
Huynen  [9]  developed  a  “ Phenomenological  Approach  to  Radar  Polarimetry” ,  which  had  a  subtle  impact 
on  the  steady  advancement  of  polarimetry  [13,  14,  15]  as  well  as  ellipsometry’  by  developing  the  “ orthogonal 
(group  theoretic)  target  scattering  matrix  decomposition ”  [28,  29,  30]  and  by  extending  the  characteristic 
optimal  polarization  state  concept  of  Kennaugh  [31,  4,  5],  which  lead  to  the  renaming  of  the  spinorial 
polarization  fork  concept  to  the  so  called  ‘ Huynen  Polarization  Fork’  in  ‘ Radar  Polarimetry  ’[31].  Here,  we 
emphasize  that  for  treating  the  general  bistatic  (asymmetric)  scattering  matrix  case,  a  more  general 
formulation  of  fundamental  Ellipsometry  and  Polarimetry  in  terms  of  a  spinorial  group-theoretic  approach  is 
strictly  required,  which  was  first  explored  by  Kennaugh  but  not  further  pursued  by  him  due  to  the  lack  of 
pertinent  mathematical  formulations  [32,  33]. 

In  ellipsometry,  the  Jones  and  Mueller  matrix  decompositions  rely  on  a  product  decomposition  of  relevant 
optical  measurement/transformation  quantities  such  as  diattenuation,  retardence,  depolarization, 
birefringence,  etc.,  [34,  35,  23,  28,  29]  measured  in  a  ‘ chain  matrix  arrangement,  i.e.,  multiplicatively 
placing  one  optical  decomposition  device  after  the  other' .  In  polarimetry,  the  Sinclair,  the  Kennaugh,  as 
well  as  the  covariance  matrix  decompositions  [29]  are  based  on  a  group-theoretic  series  expansion  in  terms 
of  the  principal  orthogonal  radar  calibration  targets  such  as  the  sphere  or  flat  plate,  the  linear  dipole  and/or 
circular  helical  scatterers,  the  dihedral  and  trihedral  comer  reflectors,  and  so  on  -  -  observed  in  a  linearly 
superimposed  aggregate  measurement  arrangement  [36,  37];  leading  to  various  canonical  target  feature 
mappings  [38]  and  sorting  as  well  as  scatter-characteristic  decomposition  theories  [39,  27,  40].  In  addition, 
polarization-dependent  speckle  and  noise  reduction  play  an  important  role  in  both  ellipsometry  and 
polarimetry,  which  in  radar  polarimetry  were  first  pursued  with  rigor  by  J-S.  Lee  [41,  42,  43,  44].  The 
implementation  of  all  of  these  novel  methods  will  fail  unless  one  is  given  fully  calibrated  scattering  matrix 
information,  which  applies  to  each  element  of  the  Jones  and  Sinclair  matrices. 
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It  is  here  noted  that  it  has  become  common  usage  to  replace  “ellipsometry”  by  “optical  polarimetry”  and 
expand  “polarimetry”  to  “radar  polarimetry”  in  order  to  avoid  confusion  [45,  18],  a  nomenclature  adopted  in 
the  remainder  of  this  paper. 

Very  remarkable  improvements  beyond  classical  “non-polarimetric”  radar  target  detection,  recognition  and 
discrimination,  and  identification  were  made  especially  with  the  introduction  of  the  covariance  matrix 
optimization  procedures  of  Tragi  [46],  Novak  et  al.  [47  -  51],  Liineburg  [52  -  55],  Cloude  [56],  and  of  Cloude 
and  Pottier  [27].  Special  attention  must  be  placed  on  the  ‘Cloude-Pottier  Polarimetric  Entropy  H , 
Anisotropy  A ,  Feature- Angle  (a  )  parametric  decomposition’  [57]  because  it  allows  for  unsupervised  target 
feature  interpretation  [57,  58].  Using  the  various  fully  polarimetric  (scattering  matrix)  target  feature 
syntheses  [59],  polarization  contrast  optimization,  [60,  61]  and  polarimetric  entropy/anisotropy  classifiers, 
very  considerable  progress  was  made  in  interpreting  and  analyzing  POL- SAR  image  features  [62,  57,  63,  64, 
65,  66].  This  includes  the  reconstruction  of  ‘Digital  Elevation  Maps  (DEMs)  ’  directly  from  POL-SAR 
Covariance-Matrix  Image  Data  Takes  ’  [67  -  69]  next  to  the  familiar  method  of  DEM  reconstruction  from 
IN-SAR  Image  data  takes  [70,  71,  72].  In  all  of  these  techniques  well  calibrated  scattering  matrix  data  takes 
are  becoming  an  essential  pre-requisite  without  which  little  can  be  achieved  [18,  19,  45,  73].  In  most  cases 
the  ‘multi-look-compressed  SAR  Image  data  take  MLC-  formatting’  suffices  also  for  completely  polarized 
SAR  image  algorithm  implementation  [74].  However,  in  the  sub-aperture  polarimetric  studies,  in 
‘ Polarimetric  SAR  Image  Data  Take  Calibration’,  and  in  ‘ POL-IN-SAR  Imaging’,  the  ‘SLC  (Single  Look 
Complex)  SAR  Image  Data  Take  Formatting’  becomes  an  absolute  must  [19,  1].  Of  course,  for  SLC- 
formatted  Image  data,  in  particular,  various  speckle-filtering  methods  must  be  applied  always. 
Implementation  of  the  ‘'Lee  Filter’  -  explored  first  by  Jong-Sen  Lee  -  for  speckle  reduction  in  polarimetric 
SAR  image  reconstruction,  and  of  the  ‘  Polarimetric  Lee-Wishart  distribution  ’  for  improving  image  feature 
characterization  have  further  contributed  toward  enhancing  the  interpretation  and  display  of  high  quality 
SAR  Imagery  [41  -  44,  75]. 


2.  The  Electromagnetic  Vector  Wave  and  Polarization  Descriptors 

The  fundamental  relations  of  radar  polarimetry  are  obtained  directly  from  Maxwell’s  equations  [86,  34], 
where  for  the  source-free  isotropic,  homogeneous,  free  space  propagation  space,  and  assuming  IEEE 
standard  [102]  time-dependence  exp(+y<ztf) ,  the  electric  E  and  magnetic  H  fields  satisfy  with  //  being  the 
free  space  permeability  and  £  the  free  space  permittivity 


VxE(r)  =  -y<z>//H(r), 
which  for  the  time-invariant  case,  result  in 

(V  +  £2)E  =  0,  E(r)  -  E0  CXP(~^T) ; 


VxEI(r)  =  jcos  E(r) 


H(r)  =  //( 


exp (-jkr) 


(2.1) 


(2.2) 


for  an  outgoing  spherical  wave  with  propagation  constant  k  =  a>  (s  //)'  and  c  =  ( £  // )  '  being  the  free 
space  velocity  of  electromagnetic  waves 

No  further  details  are  presented  here,  and  we  refer  to  Stratton  [86],  Born  and  Wolf  [34]  and  Mott  [76]  for  full 
presentations. 

2.1  Polarization  Vector  and  Complex  Polarization  Ratio 

With  the  use  of  the  standard  spherical  coordinate  system  [t~,6,(j)  ;  u;.,ue,u^j  with r,6,(f)  denoting  the 
radial,  polar,  azimuthal  coordinates,  and  U(.,U g,  U(,  the  corresponding  unit  vectors,  respectively;  the  outward 
travelling  wave  is  expressed  as 
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E  =  u,  Eg+u  E 


U  =  ueHe+u,H, 


ExH 


uf  |E|2 

2Zn 


f  V/2 

M o 

\£o  J 


=  120 n  [Q]  (2.3) 


with  P  denoting  the  Poynting  power  density  vector,  and  Z0  being  the  intrinsic  impedance  of  the  medium 

(here  vacuum).  Far  from  the  antenna  in  the  far  field  region  [86,  76],  the  radial  waves  of  (2.2)  take  on  plane 
wave  characteristics,  and  assuming  the  wave  to  travel  in  positive  z-direction  of  a  right-handed  Cartesian 
coordinate  system  (x,y,z)  ,  the  electric  field  E  ,  denoting  the  polarization  vector,  may  be  rewritten  as 


E  =  u*  Ex  +  u  v:  Ey  =|  Ex  |  cxp(  /V/y  ){at  +  uy 


exP 


(2.4) 


with  |  Ex  |,  |  E  |  being  the  amplitudes,  </)x,  <f>v  the  phases,  </>  =  </>v  - (j)x  the  relative  phase;  |  Ext Ev  |=  tana 
with  (f>x ,  <f)v ,  a  and  (j)  defining  the  Deschamps  parameters  [6,  103].  Using  these  definitions,  the 
‘ normalized  complex  polarization  vector  p  ’  and  the  ‘  complex  polarization  ratio  p  ’  can  be  defined  as 


E  u,  *  u,  I  '-  , 

|E|  |E| 


u.+ 


with  |E|~  =  E-E*  =Ex 2  +  E^  and  |E|  =  E  defines  the  wave  amplitude,  and  p  is  given  by 


(2.5) 


P  = 


<f>  =  0y  ~  </>X 


(2.6) 


2.2  The  Polarization  Ellipse  and  its  Parameters 

The  tip  of  the  real  time-varying  vector  E ,  orp,  traces  an  ellipse  for  general  phase  difference  0,  where  we 

distinguish  between  right-handed  (clockwise)  and  left-handed  (counter-clockwise)  when  viewed  by  the 
observer  in  direction  of  the  travelling  wave  [76,  19],  as  shown  in  Fig.  2.1  for  the  commonly  used  horizontal 
H  (by  replacing  x)  and  vertical  V  (by  replacing  y)  polarization  states. 

There  exist  unique  relations  between  the  alternate  representations,  as  defined  in  Fig.  2.1  and  Fig.  2.2 
with  the  definition  of  the  orientation  I//  and  ellipticity  %  angles  expressed,  respectively,  as 


\E 

a  =  \p  =  -p— 


0  <  a  <  nil 


and 


tan  2y/  =  tan(2a)  cos  (j) 


-71  !2<y/<+7t /2  (2.7) 


tan  x  -  ±minor  axis/major  axis,  sin  2 %  =  sin  2a  sin  (j)  ,  -  n  /  4  <  x  ^  n  /  4  (2.8) 

where  the  +  and  —  signs  are  for  left-  and  right-handed  polarizations  respectively. 

For  a  pair  of  orthogonal  polarizations  Pj  and  p2  =  p. 


Pi'P2*=°  A=Ai  = -1/A  » 


^l=^2+y  Zl=~Zl 


(2.9) 


In  addition,  the  following  useful  transformation  relations  exist: 
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cos2 ysm.2y/  +  /'sin2 y  ,  . 

P  = - - - -  =  tan  «  exp(y^) 

1  +  cos2jcos2^- 

where  (a,  </>)  and  (i//,  %)  are  related  by  the  following  equations: 
cos  2 a  =  cos  2^  cos  2%,  tan  <j)  =  tan  2%  /  sin  2^ 
and  inversely 


(2.10) 


(2.11) 


1 

y/  =  —  arctan 
2 


A2Re{p}A 

1  * 

1  ~PP 


+  7t  ...mod(;r)  j  =  ^-arcsin 


/2Im{p}^ 

1  * 

1  ~PP 


(2.12) 


ROTATION  SENSE:  LOOKING  INTO  THE  DIRECTION  OF  THE  WAVE  PROPAGATION 


LEFT  HANDED 

I 

ELLIPTTCIT/  ANGLE  :  x>0 


n  71 

- <  x  <  — 


RIGHT  HANDED 

I 

ELLIPTICITY  ANGLE  :  X<  0 


4  4  ^ 

(a)  Rotation  Sense  (Courtesy  of  Prof.  E.  Pottier) 


(b)  Orientation  I//  and  Ellipticity  %  Angles. 


(c)  Electric  Field  Vector. 


Fig.  2.1  Polarization  Ellipse. 
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Fig.  2.2  Polarization  Ellipse  Relations  (Courtesy  of  Prof.  E.  Pottier) 

Another  useful  formulation  of  the  polarization  vector  p  was  introduced  by  Eluynen  in  terms  of  the 
parametric  formulation  [9,  104],  derived  from  group-theoretic  considerations  based  on  the  Pauli  SU(2) 
matrix  set  y/p  {[cr]  ,  /  =  0,1, 2, 3j  as  further  pursued  by  Pottier  [105],  where  according  to  (2.10)  and 
(2.1 1),  for  y/  =  0 ,  and  then  rotating  this  ellipse  by  y/ . 


P(|E| =|  E  |  exp(j^) 


cos^  -sin^c 
sin^z  cos^ 


cos  x 
~j  sin  X 


(2.13) 


which  will  be  utilized  later  on;  and  y/p  {[cr  ]  ,  /  =  0,1, 2, 3  j  is  defined  in  terms  of  the 
classical  unitary  Pauli  matrices  [cr  ]  as 


K]  = 


i  o 
0  1 


N  = 


1  0 

0  -1 


’  K]  = 


0  1 
1  0 


’  N  = 


0  -j 
j  0 


(2.14) 


where  the  [cr/  ]  matrices  satisfy  the  unitarity  condition  as  well  as  commutation  properties  given  by 

|^{W}|=1»  ’  hlh]=K] 

satisfying  the  ordinary  matrix  product  relations. 


(2.15) 


2.3  The  Jones  Vector  and  Changes  of  Polarization  Bases 

If  instead  of  the  basis  {x  y}  or  {H  V},  we  introduce  an  alternative  presentation  {m  n}  as  a  linear  combination 
of  two  arbitrary  orthonormal  polarization  states  Em  and  E);  for  which 

E=  u.„£,„+u.,£„  (2.16) 


m  m  n  n 


and  the  standard  basis  vectors  are  in  general,  orthonormal,  i.e. 
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Vul=0,  um-ul=VuI=l  (2.17) 

with  |  denoting  the  hermitian  adjoint  operator  [21,  52,  53];  and  the  Jones  vector  Emn  may  be  defined  as 


X' 

A 1  exp  A 

=  Em 

V 

 E  exp(A) 

V 

=  |E|exp(A) 

cos  a 

A  J 

LAIexp  M_ 

V1 + pp 

L  P. 

sinaexp(A 

with  tan  a  =|  En  /  Em  |  and  (f)  =  (f)n  -  (f)m  .  This  states  that  the  Jones  vector  possesses,  in  general,  four  degrees 
of  freedom.  The  Jones  vector  descriptions  for  characteristic  polarization  states  are  provided  in  Fig.  2.3. 


[  E~ 

A 

[e.  1 

E)HII  =E(m,n)  = 

m 

A. 

E,=E(/,/)  = 

X. 

and  Eab=E(A,B)  = 

A. 

The  unique  transformation  from  the  {um  ii n }  to  the  arbitrary  {u;.  u  }  or  {ify  ii  B  }  bases  is  sought  which  is  a 
linear  transformation  in  the  two-dimensional  complex  space  so  that 

Ey  =  [U2]Emn  or  E(i,j)  =  [U2]E(m,n)  with  [U2]  [U2f  =  [I2]  (2.21) 


satisfying  wave  energy  conservation  with  [/,  ]  being  the  2x2  identity  matrix,  and  we  may  choose,  as 
shown  in  [81], 


exp(A) 

V 

-H  f.  -f,  -exP<^) 

1 

 exp  A) 

* 

-p 

^jl  +  pp 

A 

J  “  Vi  +  pT 

*-> 

A 

Jl  +  pp* 

1 

with  0'.  =  <j)j+(j)  +  7V  so  that 


1 

exp  (A) 

A  exp  (A) 

yjl  +  pp* 

/?exp(A) 

exp  (A) 

(2.22) 


(2.23) 


f 

yielding  Det{  [U2  ] }  =  exp  +  (f)'j ) }  with  <()t  +(j)j  =  0 


Since  any  monochromatic  plane  wave  can  be  expressed  as  a  linear  combination  of  two  orthonormal  linear 
polarization  states,  defining  the  reference  polarization  basis,  there  exist  an  infinite  number  of  such  bases  {i  j } 
or  [A  B}  for  which 


E  =  u„A,  +  u  A  =  A,  +  UjEj  =  u  A  +  u  A 


(2.19) 


with  corresponding  Jones  vectors  presented  in  two  alternate,  most  commonly  used  notations 


RTO-EN-SET-081 


1  -7 


Basics  of  SAR  Polarimetry  I 


ORGANIZATION 


HORIZONTAL 


H  = 


1 
0 

M f  =(/ 
X  =0 


VERTICAL 


v  = 


v= 


0 

1 

K 


2 

X  =0 


LINEAR 


ORTHOGONAL  LINEAR 


-sin  6 
cos  0 
„  7Z 

y  =  6  +  - 

2 

1  =0 


LEFT  CIRCULAR  RIGHT  CIRCULAR 


ELLIPTICAL 


ORTHOGONAL  ELLIPTICAL 


Fig.  2.3  Jones  Vector  Descriptions  for  Characteristic  Polarization  States  with  direction  of 
propagation  out  of  the  page  (Courtesy  of  Prof.  E.  Pottier) 


Since  [U2\  is  a  special  unitary  2x2  complex  matrix  with  unit  determinant,  implying  that  (i)  the  amplitude 

of  the  wave  remains  independent  of  the  change  of  the  polarization  basis,  and  that  (ii)  the  phase  of  the 
(absolute)  wave  may  be  consistently  defined  as  the  polarization  basis  is  changed,  we  finally  obtain, 


[U2]  = 


1 

1  -p 

exP  CM) 

0 

sj\  +  pp 

_p  1 

0 

exP  CM)_ 

(2.24) 


possessing  three  degrees  of  freedom  similar  to  the  normalized  Jones  vector  formulation,  but  in  most  cases  the 
phase  reference  is  taken  as  =  0  which  may  not  be  so  in  polarimetric  interferometry  [96].  For  further 
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details  on  the  group-theoretic  representations  of  the  proper  transformation  relations  see  the  formulations 
derived  by  Pottier  in  [106]. 

2.4  Complex  Polarization  Ratio  in  Different  Polarization  Bases 

Any  wave  can  be  resolved  into  two  orthogonal  components  (linearly,  circularly,  or  elliptically  polarized)  in 
the  plane  transverse  to  the  direction  of  propagation.  For  an  arbitrary  polarization  basis  {A  B}  with  unit 

vectors  a  and  b  ,  one  may  define  the  polarization  state 

YL(AB)  =  Ea  a+  Eb  b  (2.25) 


where  the  two  components  EA  and  EB  are  complex  numbers.  The  polarization  ratio  pAB  in  an  arbitrary 
basis  {A  B}  is  also  a  complex  number,  and  it  may  be  defined  as 


Pab 


Ej, 

Ea 


exp {j \(j)B  ~(pA)}  =  | PAB | exp {j(j)AB } 


(2.26) 


where  \pAB\  is  the  ratio  of  magnitude  of  two  orthogonal  components  of  the  field  EA  |  and  EB  and  (f)AH  is 

the  phase  difference  between  EA  and  EB .  The  complex  polarization  ratio  pAB  depends  on  the  polarization 
basis  {A  B}  and  can  be  used  to  specify  the  polarization  of  an  electromagnetic  wave 


E  {AB)  = 


£jexp{y^} 


1 

Pab 


\EA\^v{j<PA\ 


1 1  ebe; 
eae\ 


lebe; 


eaea 


1 

Pab 


(2.27) 


|E|exp  {j<f>A} 


1 

1 

+  Pab  Pab 

-Pab  _ 

where  |E|  =  ^E 4E*4  +  EBE*B  is  the  amplitude  of  the  wave  E (AB) .  If  we  choose  |E|  =  1  and  disregard  the 
absolute  phase  <f>A ,  the  above  representation  becomes 


E  (AB)  = 


1 

1 

Pab  Pab 

.Pab. 

(2.28) 


This  representation  of  the  polarization  state  using  the  polarization  ratio  pAB  is  very  useful.  For  example,  if 

we  want  to  represent  a  left-handed  circular  (LF1C)  polarization  state  and  a  right-handed  circular  (RF1C) 
polarization  state  in  a  linear  basis  {FI  V}  using  the  polarization  ratio.  For  a  left-handed  circular  (LFIC) 

polarization,  | EH  |  =  Ev  | ,  (j)HV  =  (f)v  -  (f)u  =  ,  and  according  to  (2.26),  the  polarization  ratio  pHV  is  j  . 

Using  (2.28)  with  pHV  =  j ,  we  obtain  for  the  left-handed  circular  (LFIC)  polarization 


E<HV)-T2 


1 

j 


(2.29) 


Similarly,  the  polarization  ratio  pHV  of  a  right-handed  circular  (RF1C)  polarization  state  in  a  linear  basis  {FI 
V}  is  -j  because  the  relative  phase  (f)HV  =  —  zr/ 2  ,  and  its  representation  is 
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E  {HV)  = 


(2.30) 


The  complex  polarization  ratio  p  is  important  in  radar  polarimetry.  However,  the  value  of  the  polarization 
ratio  p  defined  in  a  certain  polarization  basis  is  different  from  that  defined  in  the  other  polarization  basis 
even  if  the  physical  polarization  state  is  the  same. 

2.4.1  Complex  Polarization  Ratio  in  the  Linear  Basis  {H  V} 

In  the  linear  {HV}  basis  with  unit  vectors  h  and  v ,  a  polarization  state  may  be  expressed  as: 

E(HV)  =  Eh  h+  Ev  v  (2.31) 


The  polarization  ratio  pHV  ,  according  to  (2.6),  can  be  described  as: 


Phv  =' 


-‘H 


exP  ( j&Hv )  —  tan  ocHV  exp  ( j(f)HV ) ,  (f>nv  —  (j>v  (f>H 


(2.32) 


where  the  angle  aHV  is  defined  in  Fig.  2.1c,  only  in  the  {H  V}  basis  and 

\eh\  =  4Eh+ey  cos aHV 
| Ev  |  =  ^E2h  +  E2  sin  aHV 


(2.33) 


Also,  for  a  single  monochromatic,  uniform  TEM  (transverse  electromagnetic)  traveling  plane  wave  in  the 
positive  z  direction,  the  real  instantaneous  electric  field  is  written  as 


(2.34) 


£x(z,t) 

\E  x\cos(a)t  -  kz  +  </>x) 

z(z,t)  = 

£y(Z,t) 

= 

\Ey  \  cos(cot  -  kz  +  (f)Y) 

_sz{z,t)_ 

0 

In  a  cartesian  coordinate  system,  the  +x  -axis  is  commonly  chosen  as  the  horizontal  basis  (H)  and  the  +y  - 
axis  as  the  vertical  basis  (V)  Substituting  (2.33)  into  (2.34),  we  find 


£(z,t)  = 


■sJe^  +  Ey  cos  aHV  cos  {cot  - kz  +  cj)H ) 
,JeI  +  E2  sin  aHV  co s(cot  -kz  +  <j)v  ) 

\IEh+Ev  exP 


cosat 


'HV 

sin  aHV  exp (j</>) 


exp{j(cot-kz  +  </>H)} 


(2.35) 


where  (f)  =  (f)v  —  (j)H  is  the  relative  phase.  The  expression  in  the  square  bracket  is  a  spinor  [32]  which  is 
independent  of  the  time-space  dependence  of  the  traveling  wave.  The  spinor  parameters  ( a ,  if))  are  easy  to 
be  located  on  the  Poincare  sphere  and  can  be  used  to  represent  the  polarization  state  of  a  plane  wave.  In  Fig. 
2.4c,  the  polarization  state,  described  by  the  point  PE  on  the  Poincare  sphere,  can  be  expressed  in  terms  of 

these  two  angles,  where  2 aHV  is  the  angle  subtended  by  the  great  circle  drawn  from  the  point  PE  on  the 
equator  measured  from  H  toward  V;  and  (j)HV  is  the  angle  between  the  great  circle  and  the  equator. 

From  equations,  (2.7)  and  (2.8)  for  the  {H  V}  basis  we  have 
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sin  2%  =  sin  2 aHV  sin  (j)HV 
tan  2  y  =  tan(2a//K  )  cos  fiHV 


(2.36) 


which  describes  the  ellipticity  angle  %  and  the  tilt  or  orientation  angle  t//  in  terms  of  the  variables  aHV  and 
(f)HV  .  Also,  from  (2.1 1)  for  the  {H  V}  basis  an  inverse  pair  that  describes  the  aHV  and  <f)HV  in  terms  of  X 
and  y/  is  given  in  (2.37) 
cos  2  aHV  =  cos  2y/  cos  2  x 

tan2j  (2.37) 


tan  </>hv  = 


sin  2^ 


It  is  convenient  to  describe  the  polarization  state  by  either  of  the  two  set  of  angles  (aHV,<j)HV)  or  (x, tj/) 
which  describe  a  point  on  the  Poincare  sphere.  The  complex  polarization  ratio  pHV  can  be  used  to  specify 
the  polarization  of  an  electromagnetic  wave  expressed  in  the  {H  V}  basis.  Some  common  polarization  states 
expressed  in  terms  of  (^,1//)  ,  p  ,  and  the  nonnalized  Jones  vector  E  are  listed  in  Table  2.1  at  the  end  of 
this  section. 

2,4.2  Complex  Polarization  Ratio  in  the  Circular  Basis  IL  R} 

In  the  circular  basis  {L  R},  we  have  two  unit  vectors  L  (left-handed  circular)  and  R  (right-  handed 
circular).  Any  polarization  of  a  plane  wave  can  be  expressed  by 

E (LR)  =  El  L+  Er  R  (2.38) 


A  unit  amplitude  left-handed  circular  polarization  has  only  the  L  component  in  the  circular  basis  {L  R}.  It 
can  be  expressed  by 


E(ZR)  =  1*L+0*R 


1 

0 


(2.39) 


The  above  representation  of  a  unit  (LHC)  polarization  in  the  circular  basis  {L  R}  is  different  from  that  in  the 
linear  basis  {H  V}  of  (2.29).  Similarly,  a  unit  amplitude  right-handed  circular  polarization  has  only  the  R 
component  in  the  circular  basis  {L  R} 


E(LR)  =  0  *  L+ 1  *  R  = 


(2.40) 


which  is  different  from  that  in  the  linear  {H  V}  basis. 


The  polarization  ratio  pLR ,  according  to  (2.26)  is 


Plr 


E L 
El 


exp  \Mr  ~  (/>,  )}  =  |  pLR  |  exp  {j(f)IK }  =  tan  aLR  exp  {j</>LR } 


(2.41) 


where  \pLR\  is  the  ratio  of  magnitudes  of  the  two  orthogonal  components  \El  |  and  \Er  | ,  and  (f)LR  the  phase 
difference.  The  angles  aLR  and  (f)LR  are  also  easy  to  be  found  on  the  Poincare  sphere  (see  Fig.  2.6)  like  the 
angles  aHV  and  (f>HV  .  Some  common  polarization  states  in  terms  of  pLR ,  are  listed  in  Table  2.1. 


2.4.3  Complex  Polarization  Ratio  in  the  Linear  Basis  145°  135°} 

In  the  linear  {45°  135°}  basis  with  unit  vectors  45°  and  135° ,  a  polarization  state  may  be  expressed  as 
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E(45°135°)  =  E 4s„  45°  +  Em.  135°  (2.42) 

where  E^.  and  E  are  the  45°  component  and  the  135°  component,  respectively.  The  polarization  ratio 
according  to  (2.26)  is 


P. 


45°135° 


exp  {7X^35.  -  t,  )}  =  K°135°  |  ^  W 45=135"  >  =  ^  «45°135  ^  ^45=135°  >  (2‘43) 


where 


P. 


45°135” 


is  the  ratio  of  magnitudes  of  the  two  orthogonal  components 


E  . 

and 

E  . 

135° 

45° 

and  (h 

1  L 


45135° 


the  phase  difference.  The  angles  a 
Fig.  2.6) 


and  <f>45  are  also  easy  to  be  found  on  the  Poincare  sphere  (see 


TABLE  2.1 

POLARIZATION  STATES  IN  TERMS  OF  (z, '//) ,  POLARIZATION  RATIO  p  AND  NORMALIZED 

JONES  VECTOR  E 


POLARIZATION 

X 

¥ 

{H  V}  basis 

{45°  135°}  basis 

{L  R}  basis 

Phv  E 

p  E 

'45  135 

Plr  E 

Linear 

Horizontal 

0 

0 

0 

T 

0 

1  j2-r 

1  S 

T 

l 

Linear 

Vertical 

0 

n 

~2 

00 

"0" 

1 

1 

T 

l 

1  1 

1  1 

45°  Linear 

0 

n 

~4 

1  Ti 

T 

l 

0 

~l~ 

0 

lfl  ~f 
7  2 1 1  j 

135°  Linear 

0 

7C 

~4 

1 ~k 

-l 

l 

00 

"0" 

l 

.  i  r-i  -f 
7  2L-!  j. 

Left-handed 

Circular 

n 

~4 

1 

7  n 

T 

j  

1  r 1  f 

7  2  L-1  J. 

0 

~r 

0 

Right-handed 

Circular 

71 

~4 

i  rr 
~]  A-J. 

.  i r i  -f 

7  2L-1  -L 

oo 

"0" 

1 

2.5  The  Stokes  Parameters 

So  far,  we  have  seen  completely  polarized  waves  for  which \Ea  | ,  \Eb  | ,  and  (f)AH  are  constants  or  at  least 
slowly  varying  functions  of  time.  If  we  need  to  deal  with  partial  polarization,  it  is  convenient  to  use  the 
Stokes  parameters  q0,  qi ,  q2  and  q3  introduced  by  Stokes  in  1852  [107]  for  describing  partially  polarized 
waves  by  observable  power  terms  and  not  by  amplitudes  (and  phases). 

2.5.1  The  Stokes  vector  for  the  completely  polarized  wave 

For  a  monochromatic  wave,  in  the  linear  {H  V}  basis,  the  four  Stokes  parameters  are 
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K 

f+Kf 

_  i2 

\Eh\ 

2  | 

q2  =2\EH\\Ev\cos(j)HV 
<?3  =  2\Eh  | \Ey\  sin  (/)HV 


(2.44) 


For  a  completely  polarized  wave,  there  are  only  three  independent  parameters,  which  are  related  as  follows 

ql  =  q]  +  q\ + q]  (2-45) 


The  Stokes  parameters  are  sufficient  to  characterize  the  magnitude  and  the  relative  phase,  and  hence  the 
polarization  of  a  wave.  The  Stokes  parameter  q0  is  always  equal  to  the  total  power  (density)  of  the  wave;  ql 

is  equal  to  the  power  in  the  linear  horizontal  or  vertical  polarized  components;  q2  is  equal  to  the  power  in 
the  linearly  polarized  components  at  tilt  angles  y/  =  45°  or  135° ;  and  q2  is  equal  to  the  power  in  the  left- 
handed  and  right-handed  circular  polarized  components.  If  any  of  the  parameters  q0,  qx ,  q2  or  q2  has  a 
non-zero  value,  it  indicates  the  presence  of  a  polarized  component  in  the  plane  wave.  The  Stokes  parameters 
are  also  related  to  the  geometric  parameters  A  ,  %  ,  and  y/  of  the  polarization  ellipse 


<lo 

'  \eh\2+\ev\2  " 

<h 

EH\2~\Evf 

A2  cos  2 y/  cos  2% 

q2 

2  \Eh  |  \EV  cos  (fifty 

A2  cos  2^/ sin  2j 

_<h_ 

i 

c n 

<N 
_ i 

A2  sin2j 

which  for  the  normalized  case  q20  =  e2  =  e2H  +  e2  =  1  and 


9o 

\Eh\2+\Ev\2 

2  2 

e  +  e 

^ H  ^  c V 

2  1 
e 

<h 

1 

eh\2-\ev\2 

2  2 
eH  eV 

e 2  cos  2y/  cos  2% 

q2 

X/|  E„  f  +  |  Ey  I’ 

2R  e{E*HEv} 

2eHev  cos  (f) 

e 2  cos2^sin2_^ 

_<h_ 

21m  {E*hEv} 

leHev  sin  (f> 

1 

N 

<N 

.S 

X 

(N 

_ 1 

(2.46) 


(2.47) 


2.5.2  The  Stokes  vector  for  the  partially  polarized  wave 

The  Stokes  parameter  presentation  [34]  possesses  two  main  advantages  in  that  all  of  the  four  parameters  are 
measured  as  intensities,  a  crucial  fact  in  optical  polarimetry,  and  the  ability  to  present  partially  polarized 
waves  in  tenns  of  the  2x2  complex  hermitian  positive  semi-definite  wave  coherency  matrix  [./]  also  called 
the  Wolfs  coherence  matrix  [34],  defined  as: 


[J]  =  (EEt 


\e„e’„)  ( E„E ;) 

(£„£;)  (EyE-y) 


J HH 

^ HV 

1 

5 

1 _ 

qQ  +  <h  q2+j<h 
q2-j<h  q0-<h 


(2.48) 


where 


lim 

T— >oo 


indicating  temporal  or  ensemble  averaging  assuming  stationarity  of  the 


wave.  We  can  associate  the  Stokes  vector  q  with  the  coherency  matrix  [  J] 
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4o  -  \Eh\  +  \Ev  |  -  (EhEh  )  +  (EvEv  )  -  Jhh  +  Jyv 

<h  =  \EH  f  - \EV  |2  =  {EHEH  )  -  {EVEV  )  =  JHH-Jw 

'  '  '  '  (2.49) 

q2  =  2 1£„  1 |£(,  |  cos  (f>HV  =  )  =  JHV  +  JVH 

<h  =2\Eh\ \EV  I  sin  4k  =  j  {EhK  )  -  j  ( EVEH  )  =  JJHV  -  JJVH 

and  since  [  J]  is  positive  semidefmite  matrix 

Det{[J~\}>  0  or  ql>  ql  +  q;  +  q]  (2.50) 

the  diagonal  elements  presenting  the  intensities,  the  off-diagonal  elements  the  complex  cross-correlation 
between  EH  and  Ev,  and  the  Trace  {[J]} ,  representing  the  total  energy  of  the  wave.  For  JHV  =0  no 

correlation  between  EH  and  Ev  exists,  [  J\  is  diagonal  with  J HH  =  J vv  ,  (i.e.  the  wave  is  unpolarized  or 
completely  depolarized,  and  possesses  one  degree  of  freedom  only :  amplitude).  Whereas,  for 
Det{[J]}  =  0  we  find  that  JVHJHV  =  J HHJVv >  ar|d  the  correlation  between  En  and  Ev  is  maximum,  and 

the  wave  is  completely  polarized  in  which  case  the  wave  possesses  three  degrees  of  freedom:  amplitude, 
orientation,  and  ellipticity  of  the  polarization  ellipse.  Between  these  two  extreme  cases  lies  the  general  case 
of  partial  polarization,  where  Det\  [./]  j  >0  is  indicating  a  certain  degree  of  statistical  dependence  between 

Eh  and  Ev  which  can  be  expressed  in  terms  of  the  ‘ degree  of  coherency’  jtt  and  the  ‘degree  of 

polarization  ’  Dp  as 


Mhv  ~  \Mhv  exp(yAF )  — 


D 


p 


(  h1/2  /  2  2  2  \1/2 

t  4Det{[J]}  _  [<3i  +<h  +<h ) 

v  {Trace^J1)})2 )  q0 


(2.51) 


(2.52) 


where  /u  =  Dp  =  0  for  totally  depolarized  and  /u  =  Dp  =  1  for  fully  polarized  waves,  respectively. 
However,  under  a  change  of  polarization  basis  the  elements  of  the  wave  coherency  matrix  [  J]  depend  on  the 
choice  of  the  polarization  basis,  where  according  to  [52,  53],  [J]  transforms  through  a  unitary  similarity 
transformation  as 

u-53) 

The  fact  that  the  trace  and  the  determinant  of  a  hermitian  matrix  are  invariant  under  unitary  similarity 
transformations  means  that  both,  the  degree  of  polarization  as  well  as  the  total  wave  intensity  are  not 
affected  by  polarimetric  basis  transformations.  Also,  note  that  the  degree  of  coherence  flmn  does  depend  on 
the  polarization  basis.  Table  2.2  gives  the  Jones  vector  E  ,  Coherency  Matrix  [J] ,  and  Stokes  Vector  q  for 
special  cases  of  purely  monochromatic  wave  fields  in  specific  states  of  polarization. 
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TABLE  2.2 

JONES  VECTOR  E  ,  COHERENCY  MATRIX  [J] ,  AND  STOKES  VECTOR  q  FOR  SOME  STATES 

OF  POLARIZATION 


POLARIZATION 

{H  V}  BASIS 

E 

[J] 

q 

Linear 

Horizontal 

- 

f 

0 

"1  0" 

0  0 

T 

l 

0 

0 

Linear 

Vertical 

- 

0" 

i 

"0  0" 

0  1 

"  l 

-l 

0 

0 

45°  Linear 

1 

Vi 

T 

l 

'1  f 

1 

T 

0 

1 

0 

135°  Linear 

1 

VI 

"-1" 

l 

1 

2 

1  -1" 

-1  i 

"  1 

0 

-1 

0 

Left-handed 

Circular 

1 

72 

'f 

_j_ 

1 

2 

"i  -f 
J  i  _ 

T 

0 

0 

1 

Right-handed 

Circular 

i 

VI 

"  l 

-j. 

1 

2 

"  i  f 
_-J  K 

"  1 

0 

0 

-1 

2.6  The  Poincare  Polarization  Sphere 

The  Poincare  sphere,  shown  in  Fig.  2.4  for  the  representation  of  wave  polarization  using  the  Stokes  vector 
and  the  Deschamps  parameters  ( a ,  (/))  is  a  useful  graphical  aid  for  the  visualization  of  polarization  effects. 
There  is  one-to-one  correspondence  between  all  possible  polarization  states  and  points  on  the  Poincare 
sphere,  with  the  linear  polarizations  mapped  onto  the  equatorial  plane  (x  =  0)  with  the  ‘zenith  ’  presenting 
left-handed  circular  and  the  ‘nadir  ’  right-handed  circular  polarization  states  according  to  the  IEEE  standard 
notation  exp(+ jcot)  [102],  and  any  set  of  orthogonally  fully  polarized  polarization  states  being  mapped  into 
antipodal  points  on  the  Poincare  sphere  [108]. 
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EQUATOR 

LINEAR 


POLARISATION 


(a) 


Fig.  2.4  Poincare  Sphere  Representations  (Courtesy  of  Prof.  E.  Pottier) 


2.6.1  The  polarization  state  on  the  Poincare  sphere  for  the  {H  V}  basis 

In  the  Poincare  sphere  representation,  the  polarization  state  is  described  by  a  point  P  on  the  sphere,  where  the 
three  Cartesian  coordinate  components  arcr/, ,  q2 ,  and  q,  according  to  (2.46).  So,  for  any  state  of  a 

completely  polarized  wave,  there  corresponds  one  point  P(r/, ,  q^,  q3)  on  the  sphere  of  radius  q0,  and  vice 
versa.  In  Fig.  2.5,  we  can  see  that  the  longitude  and  latitude  of  the  point  P  are  related  to  the  geometric 
parameter  of  the  polarization  ellipse  and  they  are  2y/  and  2%  respectively. 
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Fig.  2.5  The  Poincare  sphere  and  the  parameters  aHV  and  (f)HV 

In  addition,  the  point  P  on  the  Poincare  sphere  can  also  be  represented  by  the  angles  aHV  and  (f)uv  .  From 
(2.37)  and  (2.46)  we  find  that 

—  =  cos2y/cos2%  =  cos2aHV  (2.54) 

<7o 

where  cos  2 aHV  is  the  direction  cosine  of  the  Stokes  vector  q  with  respect  to  the  X-axis,  i.e.,  the  angle 
2 aHV  is  the  angle  between  q  and  the  X-axis.  The  angle  (f)HV  is  the  angle  between  the  equator  and  the  great 
circle  with  basis  diameter  HV  through  the  point  P,  and  it  is  equal  to  the  angle  between  the  XOY  plane  and 
the  XOP  plane.  Drawing  a  projecting  line  from  point  P  to  the  YOZ  plane,  the  intersecting  point  P'  is  on  the 
XOP  plane,  so  (j)HV  =  Z  YO  P'  ( (j)HV  =  (f)  in  Fig.  2.5).  On  the  YOZ  plane  we  find  that 

\.WL(j)HV  =  tanZ  YOP'  =  —  =  (2.55) 

q2  sin  2y/ 

which  satisfies  equations  (2.46)  and  (2.37). 

2.6.2  The  polarization  ratio  on  the  Poincare  sphere  for  different  polarization  bases 

Also,  it  can  be  shown  that  a  polarization  state  can  be  represented  in  different  polarization  bases.  Any 
polarization  basis  consists  of  two  unit  vectors  which  are  located  at  two  corresponding  antipodal  points  on  the 
Poincare  sphere.  Fig  2.6  shows  how  the  polarization  state  P  on  the  Poincare  sphere  can  be  represented  in 
three  polarization  bases,  {HV},  {45°  135°},  and 
{L  R}.  The  complex  polarization  ratios  are  given  by 
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PhV  =\pHv\^V{.j<t>Hv) 


f tan  aHV  exp {j(j)HV) 
\~  tan  am,  cxpf  /V^,) 


a  n 

0  <  ccHV  < 


7t 

—  <a 

2 


HV 


<  rc 


P 


45“135° 


P. 


45°135° 


exP(M5n35=): 


tan  a 


45°135‘ 


exP  (7^45=135- ) 


-tan  a 


45°135' 


A  71 

0<  o'  „  <  — 

45  135  ^ 

e*P(,'0  |  <«„<, 


PlR  I  PlR  I  eXP  (  J^LR  ) 


tan  aLR  exp(y^)  °  <  a'  K  <  2 
-  tan  aLR  exp  (j(/>ut )  n_<a^<n 


(2.56) 


(2.57) 


(2.58) 


where  tan  aHV  ,  tan  «,45  l35. ,  and  tan  «//(  are  the  ratios  of  the  magnitudes  of  the  corresponding  orthogonal 
components,  and  (j)HV ,  ^  ,  and  (j)LR  are  the  phase  differences  between  the  corresponding  orthogonal 

components 


Fig  2.6  The  Polarization  State  P  in  Different  Polarization  Bases 
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2.6.3  The  relationship  between  the  Stokes  vector  and  the  polarization  ratio  for  different  polarization 

bases 

First,  consider  the  polarization  ratio  pHV  defined  in  the  {H  V}  basis.  Because  cos  2 aHV  is  the  direction 
cosine  of  the  Stokes  vector  q  with  respect  to  the  X-axis,  we  find 


_  1-tan2  aHV  J-\pHV\ 
Vo  1  +  tan  aHV  1  +  \pHV 


(2.59) 


the  straight  forward  solution  for  \pHV\  is 


\Phv\  =  , 


9o~9i 

,q0+q1 


(2.60) 


from  (2.54),  we  find 


0HV=Z  YOP'  =  tan 


(2.61) 


Combining  above  two  equations  yields 

exp  J  j  tan"1 


Phv  =  Phv  exP(j0Hv)  =  . 


q0+qx 


f  \ 

<h 

K<h; 


(2.62) 


For  a  completely  polarized  wave,  we  may  obtain  the  Stokes  vector  in  terms  of  the  polarization  ratio  pHV  by 
applying 


q0  = 

= 

q2  = 

<h  = 


V^i2  +?2  +(h  =1 


H  Phv 


2  =  cos  2  aHV 


1  + 1  Phv  | 

2 \Phv | cos $ hv  _  2  tan aHV  cos (j)HV 
1  "F  |  P Hv  |  ^  1  +  1  tan  aHV\ 

2kJsin4 


i+l/vl 


=  sin(2aWK)sin^ 


sin(2«//F)cos^//F 


(2.63) 


The  sign  of  the  three  components  of  the  Stokes  vector  is  summarized  in  Table  2.3. 

Secondly,  consider  the  polarization  ratio  , o  defined  in  the  {45°  135°}  basis.  The  cos2«4.  ^  is  the 
direction  cosine  of  the  Stokes  vector  q  with  respect  to  the  Y-axis.  So  similarly,  with 


tfo  =  1 


P 


45°135° 


^  =  tan  1 


r45°135° 


'qo+q2 
( 

_^3 

l  <h) 


(2.64) 
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TABLE  2.3 

THE  SIGN  OF  THE  qx ,  q2 ,  AND  q:]  PARAMETERS  IN  THE  {H  V}  BASIS 

</>HV 


0  <  K  <  | 


n 


aHV 

<7i 

<h 

?3 

0  <  2  a„v  <  — 

HV  2 

+ 

+ 

+ 

ft  ~ 

—  <  2a  HY  <  n 

+ 

+ 

_  _  3ft 

0  <  2  aHV  <  ^ 

' 

' 

' 

3ft  . 

—  <  2  aHV  <  2ft 

+ 

.  .  ft 

0  <  2  aHV  <  — 

+ 

+ 

ft  o 

-  <  2  aHV  <  ft 

+ 

on  3 ft 

0  <  2  aHV  <  ^ 

+ 

<  2  aHV  <  2ft 

+ 

+ 

Then  the  polarization  ratio  p45  can  be  determined  by  the  Stokes  vector  q 


P. 


9o~92 


45°135° 


exp  <  j  tan  1 


^3 

V  <hj 


\q0+q2 

Also,  the  Stokes  vector  q  can  be  determined  by  the  polarization  ratio  p  as  follows: 


?o  =1 


= 


P 


45°135° 


cos  (b 

•  c 


45°135° 


1  + 


P. 


1- 


<h  = 


45°135° 

2 


=  -Sill2«45W  C0S^- 


r45°135° 


P. 


45°135° 


1  + 


P. 


45°135° 


T  =  cos  2«45.135. 


q,= 


P, 


45°135° 


sin  6 

'  L 


45°135° 


1  + 


P< 


45°135° 


=  Sin2a45"135”  Sil1^ 


45°135° 


(2.65) 


(2.66) 


Finally,  consider  the  polarization  ratio  pLR  defined  in  the  {L  R}  basis.  Similarly,  because  the  cos  2 aLR  is 
the  direction  cosine  of  the  Stokes  vector  q  with  respect  to  the  Z-axis,  the  polarization  ratio  pLR  can  be 
determined  by  the  Stokes  vector  q  as: 


Plr 


\clo~ch 

lq0  +  q, 


exp  \  j  tan  1 


f  \ 

Uz 


(2.67) 
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Inversely, 

<7o  =  1 

_  ^ \PlR  |  C0S  ^LR  _  .  -  , 

qx  —  :  —  sm  2 ocLR  cos  (j)LR 

1  +  \Plr\ 

_2\pLR\sm</>LR  _ 

q2  2  sin  2*2^  sin  4; 

1  +  4*l 

i  l  l2 

q2  —  ■  ~2  —  cos  2  aLR 

1  +  \Plr\ 


(2.68) 


TABLE  2.4 

ALTERNATE  EXPRESSIONS  FOR  NORMALIZED  STOKES  VECTOR  PRESENTATIONS  ON  THE 

POLARIZATION  SPHERE 


X,¥ 

aHV  ’  &HV 

a  =  ',d> 

45  135  ’”45  135 

aLR  ’  &LR 

<lo 

1 

1 

1 

1 

<h 

cos  cos  2  y 

cos  2  aHV 

-Sin2«45”135»  COS^5=135- 

sin  2 aLR  cos  (f>,R 

<h 

cos  2^  sin  2  y/ 

sin  2 aHV  cos  (j)HV 

COS2°W 

sin  2 aLR  sin  (j)LR 

<h 

sin2j 

sin 2 aHV  sin  (j)nv 

Sin2«45U35=  Sil1  4=135= 

cos  2  aLR 

2.6.4  The  Poincare  polarization  sphere  and  complex  polarization  ratio  plane 

Using  the  Riemann  transformation,  Poincare  introduced  the  polarization  sphere  representation  of  Fig.  2.5 
which  gives  a  relationship  between  the  polarization  ratio  p  and  its  corresponding  spherical  coordinates  on 

the  Poincare  sphere.  First  we  need  to  introduce  an  auxiliary  complex  parameter  u(p) ,  which  is  defined  by 
the  Riemann  transformation  [14]  of  the  surface  of  the  sphere  onto  the  polar  grid  as  follows, 

n(p)  =  \^  (2.69) 

l  +  JP 

in  the  {H  V}  basis  ,  pHV  =  tan aHV  exp {j(f)HV }  =  tan  aHV  (cos  (f)uv  +  j sin <f>HV),  then 
_  (1  +  tan  aHV  sin  </>HV )  -  j  tan  aHV  cos  (f>uv 
(1  -  tan  aHV  sin  (j)HV  )  +  j  tan  aHV  cos  <j)HV 
I  ,2  _  1  +  2  tan  aHV  sin  (j)HV  )  +  tan2  aHV 
1  -  2  tan  aHV  sin  (j)HV  )  +  tan2  aHV 

4>hv 

according  to  (2.36)  and  Fig.  2.4b,  the  polar  angle  &  =  7t/2-  2%  can  be  obtained  from 
I  I2  i 

u  -1 

— 3 - =  sin  2x  =  sin  (71/  2  -  0)  =  cos  0 

14 +i 

so  that 


|u|“  —  1  2tan  aHV 

4 - =  - 4—  sin  (pHV  =  sin  2 aHV  sin 

u  +1  1  +  tan  aHV 
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0  =  cos  1 


r 

u 

2-0 

V 

u 

~+1v 

(2.70) 


also,  according  to  (2.36)  and  Fig.  2.4b,  the  spherical  azimuthal  angle  ®  =  2y/  can  be  obtained  from 
Im{u}  _  2  tan  aHV  cos  (f>HV 


Re{u} 
®  =  tan  1 


1  -  tan2  a 
Itn{u} 

Re  My 


=  tan  2 1//  =  tan  ® ,  so  that  the  spherical  azimuthal  angle  ®  becomes 


HV 


(2.71) 


Fig.  2.7  Poincare  Sphere  and  the  Complex  Plane 


2.7  Wave  Decomposition  Theorems 

The  diagonalization  of  \^JtJ  J ,  under  the  unitary  similarity  transfonnation  is  equivalent  to  finding  an 
orthonormal  polarization  basis  in  which  the  coherency  matrix  is  diagonal  or 


^  mm 

^ mn 

eu 

ei\ 

0" 

* 

en 

1 

*  04 
or 

_  ^  nm 

J nn  _ 

en 

e22_ 

0 

i 

<N 

* 

e21 

* 

e22_ 

(2.72) 

where  2,  and  A2  are  the  real  non-negative  eigenvalues  of  [./]  with  2.,  >  2,  >  0 ,  and  e,  =  [en  e12]  '  and 
c2  =  [e2l  e22  ]  1  are  the  complex  orthogonal  eigenvectors  which  define  [U2  ]  and  a  polarization  basis 

{ e, ,  e2}  in  which  [./]  is  diagonal.  [./]  is  Flermitian  and  hence  normal  and  every  normal  matrix  can  be 
unitarily  diagonalized  .  Being  positive  semidefinite  the  eigenvalues  are  nonnegative. 


2.8  The  Wave  Dichotomy  of  Partially  Polarized  Waves 

The  solution  of  (2.72)  provides  two  equivalent  interpretations  of  partially  polarized  waves  [28]:  i)  separation 
into  fully  polarized  [  JJ  ,  and  into  completely  depolarized  \J,  ]  components 


1  -22 


RTO-EN-SET-081 


Basics  of  SAR  Polarimetry  I 


[j]  =  (Ay-A2py]  +  A2[l2]  (2.72) 

where  [/2]  is  the  2x2  identity  matrix  ;  ii)  non-coherency  of  two  orthogonal  completely  polarized  wave  states 
represented  by  the  eigenvectors  and  weighed  by  their  corresponding  eigenvalues  as 
[j]  =  (/Li)[j|]+/l2[j2]  =  ^i(e,  •e1r)  +  l2(e2  -e2)  (2.74) 

where  Det{[Jy\)  =  Det{[Jf\)  =  0  ;  and  if  Ay=A2  the  wave  is  totally  depolarized  (degenerate  case) 
whereas  for  Z2  =  0 ,  the  wave  is  completely  polarized.  Both  models  are  unique  in  the  sense  that  no  other 

decomposition  in  form  of  a  separation  of  two  completely  polarized  waves  or  of  a  completely  polarized  with 
noise  is  possible  for  a  given  coherency  matrix,  which  may  be  reformulated  in  terms  of  the  ‘ degree  of 
polarization  ’  Dp  as 

Dp  =  ~Al ,  0  (A,  =  A2)  and  1  (A2  =  0)  (2.75) 

A^  +  A-, 

for  a  partially  unpolarized  and  completely  polarized  wave.  The  fact  that  the  eigenvalues  \  and  A2  are 
invariant  under  polarization  basis  transformation  makes  Dp  an  important  basis-independent  parameter. 

2.9  Polarimetric  Wave  Entropy 

Alternately  to  the  degrees  of  wave  coherency  //  and  polarization  Dp ,  the  polarimetric  wave  entropy  Hm 

[28]  provides  another  measure  of  the  correlated  wave  structure  of  the  coherency  matrix  [J] ,  where  by  using 
the  logarithmic  sum  of  eigenvalues 

H»  =  i{-PM2P,}  with  P,=-i—  (2.76) 

i=  1  Ay+A2 

so  that  Px  +  P2  =  1  and  the  normalized  wave  entropy  ranges  from  0  <  H 0)  <  1  where  for  a  completely 
polarized  wave  with  A2  =  0  and  II  (0  =  0 ,  while  a  completely  randomly  polarized  wave  with  A{  =  A2 
possesses  maximum  entropy  Hw  =  1 . 

2.10  Alternate  Formulations  of  the  Polarization  Properties  of  Electromagnetic  Vector  Waves 

There  exist  several  alternate  formulations  of  the  polarization  properties  of  electromagnetic  vector  waves 
including;  (i)  the  ‘ Four-vector  Hamiltonian’  formulation  frequently  utilized  by  Zhivotovsky  [109]  and  by 
Czyz  [110],  which  may  be  useful  in  a  more  concise  description  of  partially  polarized  waves;  (ii)  the 
‘ Spinorial  formulation  ’  used  by  Bebbington  [32],  and  in  general  gravitation  theory  [111];  and  (iii)  a  pseudo- 
spinorial  fonnulation  by  Czyz  [110]  is  in  development  which  are  most  essential  tools  for  describing  the 
general  bi-static  (non-symmetric)  scattering  matrix  cases  for  both  the  coherent  (3-D  Poincare  sphere  and  the 
3-D  polarization  spheroid)  and  the  partially  polarized  (4-D  Zhivotovsky  sphere  and  4-D  spheroid)  cases 
[109].  Because  of  the  exorbitant  excessive  additional  mathematical  tools  required,  and  not  commonly 
accessible  to  engineers  and  applied  scientists,  these  formulations  are  not  presented  here  but  deserve  our 
fullest  attention  in  future  analyses. 

3.  The  Electromagnetic  Vector  Scattering  Operator  and  the 
Polarimetric  Scattering  Matrices 

The  electromagnetic  vector  wave  interrogation  with  material  media  is  described  by  the  Scattering 
Operator  [^(k^/k,)]  with  ks,k;  representing  the  wave  vectors  of  the  scattered  and  incident, 

E'(r),E'(r)  respectively,  where 

E5 (r)  =  E s0  exp(-yks  •  r)  =  esEs0  exp(-yki  •  r)  (3.1) 

is  related  to 

E'  (r)  =  E‘0  exp(-yk,  •  r)  =  eE'()  exp(-yk .  •  r)  (3.2) 
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E* (r)  =  £XP(  j*s ' r)  [S(k,  / k . )] E' (r)  (3.3) 

r 

The  scattering  operator  [.S’ ( k  s  /k;)J  is  obtained  from  rigorous  application  of  vector  scattering  and  diffraction 

theory,  to  the  specific  scattering  scenario  under  investigation  which  is  not  further  discussed  here,  but  we  refer 
to  [97]  for  a  thought-provoking  formulation  of  these  still  open  problems. 

3.1  The  Scattering  Scenario  and  the  Scattering  Coordinate  Framework 

The  scattering  operator  [.S’ ( k  v  /  k ;. ) ]  appears  as  the  output  of  the  scattering  process  for  an  arbitrary  input  EJ,  , 

which  must  carefully  be  defined  in  terms  of  the  scattering  scenario;  and,  its  proper  unique  formulation  is  of 
intrinsic  importance  to  both  optical  and  radar  polarimetry.  Whereas  in  optical  remote  sensing  mainly  the 
1 forward  scattering  through  translucent  media  ’  is  considered,  in  radar  remote  sensing  the  ‘ back  scattering 
from  distant,  opaque  open  and  closed  surfaces'  is  of  interest,  where  in  radar  target  backscattering  we  usually 
deal  with  closed  surfaces  whereas  in  SAR  imaging  one  deals  with  open  surfaces.  These  two  distinct  cases  of 
optical  versus  radar  scattering  are  treated  separately  using  two  different  reference  frames;  the  ‘ Forward 
(anti-monostatic)  Scattering  Alignment  ( FSA )'  versus  the  'Back  Bistatic  Scattering  Alignment  ( BSA )'  from 
which  the  ' M anastatic  Arrangement '  is  derived  as  shown  in  Fig.  3.1.  In  the  following,  separately  detailed  for 
both  the  FSA  and  BSA  systems  are  shown  in  Figs.  3.2  and  3.3. 


Fig.  3.1  Comparison  of  the  FSA,  BSA,  and  MSA  Coordinate  Systems 
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z 


=  sin 9S  cos cj)sx  +  sin  9S  sin  cj)s  y  +  cos 9sz 
=  cos9s  coscj)sx  +  cos9s  sin(j)sy-sin9sz 


hf  = -sincj)sx  +  coscj)sy 

Fig.  3.2  Detailed  Forward  Scattering  Alignment  (FSA) 


z 


=  -sin9s  coscj)sx  -sin9s  sin(j)sy-cos9sz 
=  cos9s  coscj)sx  +  cos9s  sincj)sy-sin9sz 


hb  =sin^sx-cos^sy 

Fig.  3.3  Detailed  Back  Scattering  Alignment  (BSA) 


3.2  The  2x2  Jones  Forward  \J\  versus  2x2  Sinclair  [A1  Back-Scattering  Matrices 

Since  we  are  dealing  here  with  radar  polarimetry,  interferometry  and  polarimetric  interferometry,  the 
‘ bistatic  BSA  reference  frame’  is  more  relevant  and  is  here  introduced  only  for  reasons  of  brevity  but  dealing 
both  with  the  bistatic  and  the  monostatic  cases  ;  where  we  refer  to  [52,  53],  [76]  and  [19]  for  a  full  treatment 
of  the  ‘anti-monostatic  FSA  reference  frame’.  Flere,  we  refer  to  the  dissertation  of  Papathanassiou  [97],  the 
textbook  of  Mott  [76],  and  meticulous  derivations  of  Liineburg  [52]  for  more  detailed  treatments  of  the 
subject  matter,  but  we  follow  here  the  derivation  presented  in  [19].  Using  the  coordinates  of  Fig.  3.1  with 
right-handed  coordinate  systems;  xx  y,  zx ,  x2  y2  z2 ,  x3  y3  z3 ;  denoting  the  transmitter,  scatterer  and 
receiver  coordinates,  respectively,  a  wave  incident  on  the  scatterer  from  the  transmitter  is  given  by  the 
transverse  components  E  and  E  in  the  right-handed  coordinate  system  xx  yx  zx  with  the  zx  axis  pointed  at 

the  target.  The  scatterer  coordinate  system  x2  y2  z2  is  right-handed  with  z2  pointing  away  from  the  scatterer 
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toward  a  receiver.  BSA  Coordinate  System  x3  y3  z3  is  right-handed  with  z3  pointing  toward  the  scatterer.  It 
would  coincide  with  the  transmitter  system  xl  yx  zl  if  the  transmitter  and  receiver  were  co-located.  The  wave 
reflected  by  the  target  to  the  receiver  may  be  described  in  either  the  transverse  components  E  and  £  or  by 
the  reversed  components  E  and  E  .  Both  conventions  are  used,  leading  to  different  matrix  formulations.  The 
incident  and  transmitted  or  reflected  (scattered)  fields  are  given  by  two-component  vectors;  therefore  the 
relationship  between  them  must  be  a  2x2  matrix.  If  the  scattered  field  is  expressed  in  x3  y3  z3  coordinates 

(BSA),  the  fields  are  related  by  the  Sinclair  matrix  [<S] ,  thus 


Fs 

x3 

1 

S  X3X[ 

Sx}yi 

~  K 

xi 

-  j  kr2 

Fs 

\_  3  J 

J\7ir2 

_  S  yix i 

S  y^i 

i - 

i _ 

e 

(3.4) 

and  if  the  scattered  field  is  in  x2  y2  z2  coordinates  (FSA),  it  is  given  by  the  product  of  the  Jones  matrix  [./  ] 
with  the  incident  field,  thus 


ps 

1 

Tx2x  , 

T  x2y, 

K 

_  El,  _ 

V 4k  r  2 

T  y2xt 

T  y2y\ 

[  K  J 

(3.5) 


In  both  equations  the  incident  fields  are  those  at  the  target,  the  received  fields  are  measured  at  the  receiver,  and 
r2  is  the  distance  from  target  to  receiver.  The  ‘Sinclair  matrix  [/S']  ’  is  mostly  used  for  back-scattering,  but  is 
readily  extended  to  the  bistatic  scattering  case.  If  the  name  scattering  matrix  is  used  without  qualification,  it 
normally  refers  to  the  Sinclair  matrix  [5] .  In  the  general  bistatic  scattering  case,  the  elements  of  the  Sinclair 

matrix  are  not  related  to  each  other,  except  through  the  physics  of  the  scatterer.  However,  if  the  receiver  and 
transmitter  are  co-located,  as  in  the  mono-static  or  back-scattering  situation,  and  if  the  medium  between  target 
and  transmitter  is  reciprocal,  mainly  the  Sinclair  matrix  [S  (/)/?)]  is  symmetric,  i.e.  SAB  =  SBA .  The  Jones 

matrix  is  used  for  the  forward  transmission  case;  and  if  the  medium  between  target  and  transmitter,  without 
Faraday  rotation,  the  Jones  matrix  is  usually  normal.  However,  it  should  be  noted  that  the  Jones  matrix  is  not  in 
general  normal,  i.e.,  in  general  the  Jones  matrix  does  not  have  orthogonal  eigenvectors.  Even  the  case  of  only 
one  eigenvector  (and  a  generalized  eigenvector)  has  been  considered  in  optics  (homogeneous  and 
inhomogeneous  Jones  matrices).  If  the  coordinate  systems  being  used  are  kept  in  mind,  the  numerical 
subscripts  can  be  dropped. 

It  is  clear  that  in  the  bistatic  case,  the  matrix  elements  for  a  target  depend  on  the  orientation  of  the  target  with 
respect  to  the  line  of  sight  from  transmitter  to  target  and  on  its  orientation  with  respect  to  the  target-receiver  line 
of  sight.  In  the  forms  (3.4)  and  (3.5)  the  matrices  are  absolute  matrices,  and  with  their  use  the  phase  of  the 
scattered  wave  can  be  related  to  the  phase  of  the  transmitted  wave,  which  is  strictly  required  in  the  case  of 
polarimetric  interferometry.  If  this  phase  relationship  is  of  no  interest,  as  in  the  case  of  mono-static  polarimetry, 
the  distinct  phase  term  can  be  neglected,  and  one  of  the  matrix  elements  can  be  taken  as  real.  The  resulting  form 
of  the  Sinclair  matrix  is  called  the  relative  scattering  matrix.  In  general  the  elements  of  the  scattering  matrix 
are  dependent  on  the  frequency  of  the  illuminating  wave  [19,  14,  15]. 

Another  target  matrix  parameter  that  should  be  familiar  to  all  who  are  interested  in  microwave  remote  sensing  is 
the  radar  cross  section  (RCS).  It  is  proportional  to  the  power  received  by  a  radar  and  is  the  area  of  an 
equivalent  target  that  intercepts  a  power  equal  to  its  area  multiplied  by  the  power  density  of  an  incident  wave  and 
re-radiates  it  equally  in  all  directions  to  yield  a  receiver  power  equal  to  that  produced  by  the  real  target.  The  radar 
cross  section  depends  on  the  polarization  of  both  transmitting  and  receiving  antennas.  Thus  the  radar  cross 
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section  may  be  specified  as  HH  (horizontal  receiving  and  transmitting  antennas),  HV  (horizontal  receiving  and 
vertical  transmitting  antennas),  etc.  When  considering  ground  reflections,  the  cross  section  is  normalized  by  the 
size  of  the  ground  patch  illuminated  by  the  wave  from  the  radar.  The  cross  section  is  not  sufficient  to  describe 

the  polarimetric  behavior  of  a  target.  In  terms  of  the  Sinclair  matrix  [.S'] ,  and  the  normalized  effective  lengths  of 
transmitting  and  receiving  antennas,  h  and  h  ,  respectively,  the  radar  cross  section  is 


(7  rt 


h,T  [S]  h, 


(3.6) 


A  polarimetrically  correct  form  of  the  radar  equation  that  specifies  received  power  in  terms  of  antenna  and 
target  parameters  is 


Wrt 


W,Gj0,4)Ae/0,4) 


(4n:  nr 2T 


hrT[S]ht 


(3.7) 


where  Wt  is  the  transmitter  power  and  subscripts  t  and  r  identify  transmitter  and  receiver,  and  its  properties  are 
defined  in  more  detail  in  Mott  [76]  and  in  [19].  The  effective  antenna  height  h (6?  ,</>),  is  defined  via  the  electric 
field  E*  (r,  6  ,  (f),  radiated  by  an  antenna  in  its  far  field,  as 

E( (r, 0,</>)=  - exp(-  jkr)h(<9  ,</>),  (3.8) 

2  A  r 

with  Z0  the  characteristic  impedance,  A,  the  wavelength,  and  /  the  antenna  current. 

3.3  Basis  Transformations  of  the  2x2  Sinclair  Scattering  Matrix  [5] 

Redefining  the  incident  and  scattering  cases  in  terms  of  the  standard  {H  V}  notation  with  H  =  x,  V  =  y 
and  with  proper  re-normalization,  we  redefine  (3.1)  as 

Vhv=[S]hvKv  or  K\HV)  =  [S(HV)]E\HV)  (3.9) 

where  the  complex  conjugation  results  from  inversion  of  the  coordinate  system  in  the  BSA  arrangement 
which  invites  a  more  rigorous  formulation  in  terms  of  directional  Sinclair  vectors  including  the  concepts  of 
time  reversal  as  treated  by  Liineburg  [52].  Using  these  Sinclair  vector  definitions  one  can  show  that  the 
transformation  from  one  orthogonal  polarisation  basis 

{H  V}  into  another  {i  j }  or  [A  B}  is  a  unitary  congruence  (unitary  consimilarity)  transformation  of  the 
original  Sinclair  scattering  matrix  [.S’]//f/  into  [.S  ],; ,  where 

[Slj  =[C/2][5]^[C/2]  r  or  [S(//)]  =  [t/2]  [S{HV)]  =  [U2]  t  (3.10) 

with  [f/2  ]  given  by  (2.23),  so  that  the  components  of  the  general  non-symmetric  scattering  matrix  for  the 
bistatic  case  in  the  new  polarization  basis,  characterized  by  a  complex  polarization  ratio  p  ,  can  be  written  as 
[81,25] 

S,  =  — ! -pS„-pSy„  +p2S„ ]  (3.11) 

1 +PP 

S„=—i-^\pSm*Sln,-pp-SrH-pS„\ 

\+pp 
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Sj,  =  km  -  pps„y  +  S,.„  -  ps„  ] 

1 +  PP 

Sjj  =  i  +  pp*  \-P  Shh  +  PSHV  ~  PSVH  +  ^w] 

There  exist  three  invariants  for  the  general  bistatic  case  (BSA)  under  the  change-of-basis  transformation  as 
given  by  (3.5): 


(i)  k4  =  Span[S]  =  SHH  \z  +  \  SHV  \z  +  \  SVH  \z  +  \  Svv  |>  Su  \z  +  \  S,  \2  +  \  S„  \ i  +  \  S, 


(3.12) 


confirms  that  the  total  power  is  conserved,  and  it  is  known  as  Kennaugh’s  span-invariant  k4  ; 

(ii)  SHV  —SVH  =Sj  —  S jj  ,  for  monostatic  case  (3.13) 

warranting  symmetry  of  the  scattering  matrix  in  any  polarization  basis  as  long  as  the  BSA  for  the  strictly 
mono-static  but  not  general  bistatic  case  is  implied; 

(iii)  Det{[S]HV}  =  Det{[S\}  or  Det { [S (HV)] }  =  Det { [5* ((/)] }  (3.14) 

due  to  the  fact  that  Det{[(/2]}  =  1  implies  determinantal  invariance. 

In  addition,  diagonalization  of  the  scattering  matrix,  also  for  the  general  bistatic  case,  can  always  be  obtained 
but  requires  mixed  basis  representations  by  using  the  ‘Singular  Value  Decomposition  Theorem  (SVD)  ’  [52, 
53]  so  that  the  diagonalized  scattering  matrix  [5),]  can  be  obtained  by  the  left  and  right  singular  vectors, 
where 


[Sd]  =  [Ql][S][Qs]  with  [SD]  = 


and  Z)ef{[0J}  =  £><*{[&]}  =1 


(3.15) 


and  and  Sy  denote  the  diagonal  eigenvalues  of  [5] ,  and  the  diagonal  elements  SA  and  Sy  can  be 

taken  as  real  nonnegative  and  are  known  as  the  singular  values  of  the  matrix  [5] .  For  the  symmetric 
scattering  matrices  in  the  mono-static  case  (MSA),  diagonalization  is  achieved  according  to  the  unitary 
consimilarity  transform  for  which 

[Qr\  =  [QlY  (3.16) 


and  above  equations  will  simplify  due  to  the  restriction  of  symmetry  5..  =  Sjt . 

3.4  The  4x4  Mueller  (Forward  Scattering)  [Ml  and  the  4x4  Kennaugh  (Back-Scattering)  \K\ 

Power  Density  Matrices 

For  the  partially  polarized  cases,  there  exists  an  alternate  formulation  of  expressing  the  scattered  wave  in 
terms  of  the  incident  wave  via  the  4x4  Mueller [M]  and  Kennaugh  [A  ]  matrices  for  the  FSA  and  BSA 
coordinate  formulations,  respectively,  where 

[q’]  =  [M][q‘]  (3.17) 
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For  the  purely  coherent  case,  [M]  can  formally  be  related  to  the  coherent  Jones  Scattering  Matrix  [T]  as 
[M]  =  [  1  1  1  - 1]  [Af~l  ( [r]  0  |Y]*  )[All  =  [A](  [r]  0  [T]*  )[A\l  (S.18) 


with  ®  symbolizing  the  standard  Kronecker  tensorial  matrix  product  relations  [112]  provided  in  (A.  1), 
Appendix  A,  and  the  4x4  expansion  matrix  [A]  is  given  by  [76]  as 


1  0 
1  0 
0  1 
0  j 


0  1 
0  -1 
1  0 

~j  0 


(3.19) 


with  the  elements  M..  of[M],  given  in  (B.l),  Appendix  B. 

Specifically  we  find  that  if  [J]  is  normal,  i.e.  [r][r]r  =  [T]r  [7]  ,  then  [M]  is  also  normal,  i.e. 
[M][M]t  =  [M]t[M]  . 

Similarly,  for  the  purely  coherent  case  [76],  [A]  can  formally  be  related  to  the  coherent  Sinclair 
matrix  [ S ]  as 

[if]  =  2[A]ri([1S]®[1S]*)[A]1  (3.20) 

where 


[A]™=±[AT  (3.21) 

and  for  a  symmetric  Sinclair  matrix  [.S']  ,  then  [A]  is  symmetric,  keeping  in  mind  the  ‘ mathematical 
formalism’  [M]  =  diag[  1  1  1  -  l][/v ] ,  but  great  care  must  be  taken  in  strictly  distinguishing  the 
physical  meaning  of  [A]  versus  [M]  in  terms  of  [S]  versus  [T]  respectively.  Thus,  if  [S]  is 
symmetric,  SHV  =  SVH  ,  then  [if]  is  symmetric,  AT  =  K ..  ;  and  the  correct  elements  for  [M] ,  [A"]  and 
the  symmetric  cases  are  presented  in  (B.l  -  B.5),  Appendix  B. 

3.5  The  2x2  Graves  Polarization  Power  Scattering  Matrix  \G 1 

Kennaugh  introduces,  next  to  the  Kennaugh  matrix  [A'] ,  another  formulation  [G] ,  for  expressing  the  power 

in  the  scattered  wave  E,s  to  the  incident  wave  E'  for  the  coherent  case  in  terms  of  the  so-called  ‘ Graves 
polarization  coherent  power  scattering  matrix  [G]  ’,  where 

PS  =~  ~  jEiT*[G]Ei  (3.22) 

8*  Z 0r~ 

so  that  in  terms  of  the  Kennaugh  elements  AT ,  defined  in  the  appendix,  for  the  mono-static  case 
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[G]=<mr[s]> 


Kn+Kn  Ku-jKu 

K\3  +  jK  14  K1 ,  —  Kn 


(3.23) 


By  using  a  single  coordinate  system  for  (x1,y1,zl)  and  (x3,  v3,z3)  for  the  monostatic  case,  as  in  Fig.  3.1, 

and  also  described  in  detail  in  [19],  it  can  be  shown  that  for  a  scatterer  ensemble  (e.g.  precipitation)  for 
which  individual  scatterers  move  slowly  compared  to  a  period  of  the  illuminating  wave,  and  quickly 
compared  to  the  time-averaging  of  the  receiver,  time-averaging  can  be  adjusted  to  find  the  decomposed 
power  scattering  matrix  ([g])  ,  as 


([G])  =  <[S(Of  [S(0]>  =  [G„]+[Gk] 


|2 

* 

HV 


ShhS 

I  sHV 


HV 

\2 


+  < 


I  ^VH  I 

9  O* 

^VH^VV 


(3.24) 


This  shows  that  the  time  averaged  ‘Graves  Power  Scattering  Matrix  ([G])  ’,  first  introduced  by  Kennaugh  [4, 

5],  can  be  used  to  divide  the  powers  that  are  received  by  linear  horizontally  and  vertically  by  polarized 
antennas,  as  discussed  in  more  detail  in  [19]  and  in  [113].  It  should  be  noted  that  a  similar  decomposition 
also  exists  for  the  Muller/Jones  matrices,  commonly  denoted  as  FSA  power  scattering  matrix 


Thh^HV 


Thh^hv 


gVH  I 

^VH^VV 


^VH^VV 


lvv\ 


(3.25) 


which  is  not  further  analyzed  here  [113]. 


3.6  Co/Cross-Polar  Backscattering  Power  Decomposition  for  the  One-Antenna  (Transceiver)  and 

the  Matched  Two-Antenna  (Quasi-Monostatic)  Cases 

Assuming  that  the  scatterer  is  placed  in  free  unbounded  space  and  that  no  polarization  state  transformation 
occurs  along  the  propagation  path  from  the  transmitter  (T)  to  the  scatterer  incidence  (S),  and  along  that  from 
the  scatterer(s)  to  the  receiver  (R),  then  the  value  of  the  terminal  voltage  of  the  receiver,  VR  ,  induced  by  an 
arbitrarily  scattered  wave  E/(  at  the  receiver,  is  defined  by  the  radar  brightness  function  VR ,  and  the 
corresponding  received  power  PR  expression 


rx=K** 


PR  V« 


(3.26) 


with  the  definition  of  the  Kennaugh  matrix  [A-]  in  terms  of  the  Sinclair  matrix  [5] ,  the  received  power  or 
radar  brightness  function  may  be  re-expressed 


p 

'  RT 


h,  [S]  e7 


1  T 

^  Qfi  [^]  Qr 


(3.27) 


where  q/f  and  qr  the  corresponding  normalized  Stokes’  vectors. 

For  the  one-antenna  (transceiver)  case  the  co-polar  channel  (c)  and  the  cross-polar  channel  (x)  powers 
become: 


P=L 

c  2 


hr  [S]E3 


2  1  T 

=-q t  [kmt 


(3.28) 
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*=  2 

with 


hr/  [S]Er  =  |  q/  [Iv]  qr 


(3.29) 


[Kc  ]  =  {[A]  1  j  ([F]  0  [Tf  )[AY 1  =  [C]  [f] 

]  =  M'1  J  (M  [T]  0  [F]*  )[A]1  =  [Cx  ]  [K] 


(3.30) 

(3.31) 


and 
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(3.32) 


For  the  Two-Antenna  Dual  Polarization  case,  in  which  one  antenna  serves  as  a  transmitter  and  the  other  as 
the  receiver,  the  optimal  received  power  Pm  for  the  ‘matched  case '  becomes  by  using  the  matching  condition 


V  =W. 


so  that 


(3.33) 


„]qr,  where  [XJ  =  fc]+[^]=[^,]h],and  [A'n]  = 


10  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 


(3.34) 


which  represent  an  essential  relationship  for  determining  the  optimal  polarization  states  from  the 
optimization  of  the  Kennaugh  matrix. 

3.7  The  Scattering  Feature  Vectors  :  The  Lexicographic  and  the  Pauli  Feature  Vectors 

Up  to  now  we  have  introduced  three  descriptions  of  the  scattering  processes  in  terms  of  the  2x2  Jones  versus 
Sinclair,  [T]  versus  [.S'] ,  the  2x2  power  scattering  matrices,  [F]  versus  [G] ,  and  the  4x4  power  density 
Muller  versus  Kennaugh  matrices,  [M]  versus  [F] .  Alternatively,  the  polarimetric  scattering  problem  can 
be  addressed  in  terms  of  a  vectorial  feature  descriptive  formulation  [114]  borrowed  from  vector  signal 
estimation  theory.  This  approach  replaces  the  2x2  scattering  matrices  [F]  versus  [.S'] ,  the  2x2  power 
scattering  matrices  [F]  versus  [G] ,  and  the  4x4  Muller  [M]  versus  Kennaugh  [F]  matrices  by  an 
equivalent  four-dimensional  complex  scattering  feature  vector  f4 ,  formally  defined  for  the  general  bi-static 
case  as 


[S]hv 


f,  =  F{[S]}  =  ^Trace{[S]ip}  =  [f0  fx  f2  /3f 


(3.35) 


where  F{  [iS] }  is  the  matrix  vectorization  operator  Trace{[S] }  is  the  sum  of  the  diagonal  elements  of  [.S] , 
and  iff  is  a  complete  set  of  2x2  complex  basis  matrices  under  a  hermitian  inner  product.  For  the 
vectorization  of  any  complete  orthonormal  basis  set  [97]  of  four  2x2  matrices  that  leave  the  (Euclidean) 
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norm  of  the  scattering  feature  vector  invariant,  can  be  used,  and  there  are  two  such  bases  favored  in  the 
polarimetric  radar  literature;  one  being  the  ‘lexicographic  basis  ,  and  the  other  ‘Pauli  spin  matrix  set  ’ 

hvf  We  note  here  that  the  distinction  between  the  lexicographic  and  Pauli-based  feature  vector 

representation  is  related  to  Principal  and  Independent  Component  Analysis  (PCA/ICA)  which  is  an 
interesting  topic  for  future  research. 

(i)  The  ‘Lexicographic  Feature  vector  f4/ ,  is  obtained  from  the  simple  lexicographic  expansion  of  [5] 
using  P'd  ,  with 


so  that  the  corresponding  feature  vector  becomes 

Ur.  =  [Srh  Shy  Svh  Svv  ]r  (3.37) 


1  0 
0  0 


0  1 
0  0 


0  0 
1  0 


0  0 
0  1 


(3.36) 


(ii)  The  Pauli  Feature  vector  Up  is  obtained  from  the  renowned  complex  Pauli  spin  matrix  basis  set 
which  in  a  properly  re-normalized  presentation  is  here  defined  as 


frd- \F 


1  0 
0  1 


v/2 


1  0 

0  -1 


v/2 


0  1 
1  0 


v/2 


0  -j 

j  0 


(3.38) 


resulting  in  the  ‘polarimetric  correlation  phase  ’  preserving  ‘Pauli  Feature  Vector 

f4 .P=[/o  f,  A  A\F^[Sm+S„  S„„-S„h  S„r+Sm  j(SH„-Sm)]T  (3.39) 


where  the  corresponding  scattering  matrix  [5*]^  is  related  to  the  f4P  =  [/0  f  f2  /3] 


by 


=  F 


fo  ~fi  fi  ~jfi 
fl  +  jf\  fo  +  f\ 


=F1 


(3.40) 


3.8  The  Unitary  Transformations  of  the  Feature  Vectors 

The  insertion  of  the  factor  2  in  (3.36)  versus  the  factor  V2  in  (3.38)  arises  from  the  ‘total  power 
invariance' ,  i.e.  keeping  the  norm  independent  from  the  choice  of  the  basis  matrices  VF  ,  so  that 

WK  f<  =  \ Span{[s]}=2rrac4-S]h]f}=  +  IVI2  +I5™!'  +|V|2)=  k4  (3.41) 

This  constraint  forces  the  transformation  from  the  lexicographic  to  the  Pauli-based  feature  vector  [52,  53, 
1 14],  or  to  any  other  desirable  one,  to  be  unitary,  where  with 

f4P=[A,]f4L  and  reversely  fn.  =  [aJ'^p  (3-42) 

we  find 
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[D4]  = 


1 

V2 


"1 

0 

0 

1 ' 

"1 

1 

0 

0  " 

1 

0 

0 

-1 

,  r  n,  1 

0 

0 

1 

-  j 

0 

1 

1 

0 

'l^ 

II 

II 

T 

■'4' 

Q 

0 

0 

1 

j 

0 

j 

-j 

0 

1 

-1 

0 

0 

(3.43) 


Furthermore,  these  special  unitary  matrices  relating  the  feature  vectors  control  the  more  general  cases  of 
transformations  related  to  the  change  of  polarization  basis.  By  employing  the  Kronecker  direct  tensorial 
product  of  matrices,  symbolized  by  0 ,  we  obtain,  the  transformation  for  the  scattering  vector  from  the 


linear  {uw,uK}  to  any  other  elliptical  polarization  basis  j u A ,u B j ,  characterized  by  the  complex 


polarization  ratio  by 


f4z (4B)  =  [U 4L ]f4i (HV)  and  f4P(AB)  =  [U4P]f4P(HV) 


(3.44) 


where  [U4L  ]  is  the  transformation  matrix  for  the  conventional  feature  vector  f4£ 

Flere  we  note  that  in  order  to  obtain  the  expression  [U4L  ]  =  [U2  ]  0  [f/,]7  ,  the  unitary  congruence  (unitary 

consimilarity)  transformation  for  the  Sinclair  scattering  matrix  in  the  reciprocal  case  was  used.  This  implies 
however  that  we  must  distinguish  between  forward  scattering  and  backscattering  (and  so  also  bistatic 
scattering);  where  for  the  reciprocal  backscatter  case  the  3-dimensional  target  feature  vectors  ought  to  be 
used.  These  features  lead  to  interesting  questions  which  however  need  more  in  depth  analyses  for  which  the 
ubiquity  of  the  Time  Reversal  operation  shows  up  again. 


[U4L]  =  [U2]®[U2f 


1 

-I  * 

l+pp 


1 

p 

p 


*  * 

-p  -  p  p 

1  - pp *  -p 

-pp  1  -p 
p  p  1 


and  \U4P  ]  is  the  homologous  transformation  matrix  for  the  Pauli-based  feature  vector  f4P 


(3.45) 


[U4P]  =  [D4][U4L][D4f 


(3.46) 


where  [U4L  \  and  \U4P]  are  special  4x4  unitary  matrices  for  which  with  [/4]  denoting  the  4x4  identity 
matrix 


[U4][U4]  =  [I4]  and  Det{[U4]}  =  1  (3.47) 

Kennaugh  matrices  and  covariance  matrices  are  based  on  completely  different  concepts  (notwithstanding 
their  formal  relationships)  and  must  be  clearly  separated  which  is  another  topic  for  future  research. 

The  main  advantage  of  using  the  scattering  feature  vector,  f4£  or  f4P,  instead  of  the  Sinclair  scattering 
matrix  [5]  and  the  Kennaugh  matrix  [K] ,  is  that  it  enables  the  introduction  of  the  covariance  matrix 
decomposition  for  partial  scatterers  of  a  dynamic  scattering  environment.  Flowever,  there  does  not  exist  a 
physical  but  only  a  strict  relationship  mathematical  between  the  two  alternate  concepts  for  treating  the 
partially  coherent  case,  which  is  established  and  needs  always  to  be  kept  in  focus  [114],  It  should  be  noted 
that  besides  the  covariance  matrices  the  so-called  (normalized)  correlation  matrices  are  often  used 
advantageously  especially  when  the  eigenvalues  of  a  covariance  matrix  have  large  variations. 

3.9  The  Polarimetric  Covariance  Matrix 

In  most  radar  applications,  the  scatterers  are  situated  in  a  dynamically  changing  environment  and  are  subject 
to  spatial  (different  view  angles  as  in  ‘SAR’)  and  temporal  variations  (different  hydro-meteoric  states  in 
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‘RAD-MET’),  if  when  illuminated  by  a  monochromatic  waves  cause  the  back-scattered  wave  to  be  partially 
polarized  with  incoherent  scattering  contributions  so  that  “([5*])  =  {[5(r, /)]} Such  scatterers,  analogous 

to  the  partially  polarized  waves  are  called  partial  scatterers  [78,  90].  Whereas  the  Stokes  vector,  the  wave 
coherency  matrix,  and  the  Kennaugh/Mueller  matrix  representations  provide  a  first  approach  of  dealing  with 
partial  scattering  descriptions,  the  unitary  matrix  derived  from  the  scattering  feature  f4  vector  provides 

another  approach  borrowed  from  decision  and  estimation  signal  theory  [115]  which  are  currently  introduced 
in  Polarimetric  SAR  and  Polarimetric-Interferometric  SAR  analyses,  and  these  need  to  be  introduced  here. 
However,  even  if  the  environment  is  dynamically  changing  one  has  to  make  assumption  concerning 
stationarity  (at  least  over  timescales  of  interest),  homogeneity  and  ergodicity.  This  can  be  analyzed  more 
precise  by  introducing  the  concept  of  space  and  time  varying  stochastic  processes. 

The  4x4  lexicographic  polarimetric  covariance  matrix  [C4i]  and  the  Pauli-based  covariance  matrix 
[C4P]  are  defined,  using  the  outer  product  ®  of  the  feature  vector  with  its  conjugate  transpose  as 


\CAL\  =  <f4i  -fJi)  and  [C4P]  =  (f4P  -fjp) 


(3.48) 


where  (...'j  indicates  temporal  or  spatial  ensemble  averaging,  assuming  homogeneity  of  the  random  medium. 
The  lexicographic  covariance  matrix  [C4]  contains  the  complete  information  in  amplitude  and  phase  variance 
and  correlation  for  all  complex  elements  of  [5]  with 


(3.49) 


and 

[C4P]  =  (f4P-0  =  <[Z>4]f4i  f tL[Dj)  =  [D4]( f4L  f :L)\D4f  =  [D4][C4L][D4f  (3.50) 

Both  the  ‘Lexicographic  Covariance  [C4/  ]  ‘  and  the  ‘ Pauli-based  Covariance  [C4/J]  /  ’  matrices  are 

hermitian  positive  semi-definite  matrices  which  implies  that  these  possess  real  non-negative  eigenvalues  and 
orthogonal  eigenvectors.  Incidentally,  those  can  be  mathematically  related  directly  to  the  Kennaugh  matrix 
[A],  which  is  not  shown  here;  however,  there  does  not  exist  a  physical  relationship  between  the  two 
presentations  which  must  always  be  kept  in  focus. 

The  transition  of  the  covariance  matrix  from  the  particular  linear  polarization  reference  basis 

{H  V}  into  another  elliptical  basis  [A  B},  using  the  change-of-basis  transformations  defined  in  (3.41  -  3.45), 

where  for 

[C4L{AB)\  =  (i4L{AB)-ftL{AB))  =  [U, 4  ](f4i  (HV)  ■  f4i  (HV))[U  4  ]+  =  [D4][C4i(^F)][D4]t  (3.51) 

and  for 

[c4AAB)]=(f4P(AB)-fi(AB))=[u4](f4Affn-fU^m^4]f=[mc4A^nm]f 
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The  lexicographic  and  Pauli-based  covariance  matrices,  [C4i]  and  [C4P] ,  contain,  in  the  most  general  case, 

as  defined  in  (3.49)  and  (3.50),  sixteen  independent  parameters,  namely  four  real  power  densities  and  six 
complex  phase  correlation  parameters. 

3.10  The  Monostatic  Reciprocal  Back-Scattering  Cases 

For  a  reciprocal  target  matrix,  in  the  mono-static  (backscattering)  case,  the  reciprocity  constrains  the  Jones 
matrix  to  be  usually  normal,  and  the  Sinclair  scattering  matrix  to  be  symmetrical,  i.e.  SHV  =  SVH  ,  which 
further  reduces  the  expressions  of  [G]  and  [A-] .  Furthermore,  the  four-dimensional  scattering  feature  vector 
f4  reduces  to  a  three-  dimensional  scattering  feature  vector  f3  such  that  following  [97] 

f 3Z.  =[0]  >  f«  = 

where  with  [/3]  being  the  unit  3x3  matrix  and  always  keeping  in  mind  that  the  transformation  between 

lexicographic  and  Pauli  ordering  is  a  direct  transformation  of  the  scattering  matrix  (and  not  only  of  the 
covariance  matrices) 


yf2Sh 


Sm  =  S„ 


(3.53) 


[0]  = 


10  0  0 

0  '  '  0 

V2  V2 

0  0  0  1 


and  [Q][Q]r  =[/3] 


(3.54) 


and  the  factor  VI  needs  to  be  retained  in  order  to  keep  the  vector  norm  consistent  with  the  span  invariance 

K . 

Similarly,  the  reduced  Pauli  feature  vector  f3p  becomes 


f3P  =[Q\  J 


4  P 


-I-  9 

HH  ^  UW 


C  _  C 

UHH  ^VV 


2  V] 


c  _  e 

^HV  UVH 


(3.55) 


The  three-dimensional  scattering  feature  vector  from  the  lexicographic  to  the  Pauli-based  matrix  basis,  and 
vice  versa,  are  related  as 

f3p=[G3]f3i  and  f3i=[03r1f3P  (3.56) 

with  [D3]  defining  a  special  3x3  unitary  matrix 


m= 


i 

vi 


“l 

0 

l  “ 

and  [D3]  '  =[ZJ3]:  =1 

"1 
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0  " 
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0 

-l 

0 
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VI 
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1  -1 

0 

(3.57) 


The  change-of-basis  transformation  for  the  reduced  scattering  vectors  in  terms  of  the  complex  polarization 
ratio  p  of  the  new  basis  is  given  by 

i,L{AB)  =  [U,L{p)JiL{HV)  and  f3P (AB)  =  [U3P (p)]f3P (HV)  (3.58) 

where 
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\U3L\  = 


1 +PP* 


1  42  p  p2 

-  42  p*  1  -  pp *  42  p 

p2  -42p  1 


(3.59) 


and 


[U3P]  =  [D3][U3L][D3] 


t 


I 


2(1  +  PP) 


2  +  p2  +  p* 

2  *2 

P  ~p 

2  ip~p) 


* 


p  -p~  2  {p-p) 

2 -{p2+p4  2 (p  +  p) 

~2  (p  +  p)  2(1  -pp) 


(3.60) 


which  are  3x3  special  unitary  matrices. 

Thus,  a  reciprocal  scatterer  is  completely  described  either  by  the  3x3  ‘Polarimetric  Covariance  Matrix 

[C,  J’ 


J  HH 


Ms„A)  ( sHHs;„ 


•J2(sms-)  2l\S 


J  HV^HH 

c  SI* 

^VV^HH 


HV 

42(svvs 


VV^HV 


42(s 
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HV^VV 


vv 


(3.61) 


or  by  the  3x3  Polarimetric  Pauli  Coherency  Matrix  [C3P] 


[c„]=< 


(3.62) 


C  -I-  ^ 

^  HH  '  ^VV 


(sHH-svv)(sHH+svvy 

2(Shv(Shh  +Svv)*) 


(SHH+SVV)(SHH - svvy )  2((shh  +sw)s*Hr ) 
2)  2((Shh  -SVV)S*HV) 

2(sHV(sHH-srvy)  -/ip  l2' 


Shh  Syr 


4  15 


HV 


where  the  relation  between  the  3x3  Pauli  coherency  matrix  [C3P]  and  the  3x3  covariance  matrix  [C3/  ]  is 
given  by 

[C3L\  =  <f3z  O  =  ([A]f3p  =  [A]<f3P-f3P>[A]+  =  [A][C3p][A]+  (3-63) 

and 

[C3i  (AB)]  =  (f3L  (AB)  ■  A (AB))  =  [U3L  ]<f3i  (HV)  •  A  (HV))[U3L  ]+  =  [U3L ][C3Z (HV)] [U3L ]+  (3.64) 


\C3P(AB)\  =  (f3p(AB)Al(AB))  =  [U3P](f3AHV)Al(HV))[U3py=[U3P][C3P(HV)][U3Fr(3.65) 


wher  e[U3L(p)][U3L(p)]+ =[I3]  and  Det{[U3L(p)]}  =  1 


(3.66) 


and 


HA!2  =  liyf  =  \sPan{[S]}  =  ^Trace{[S][Sf}  =  ±  {|Sj2  +2|5j2  *|y2 


=  k3  (3.67) 
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3.11  Co/Cross-polar  Power  Density  and  Phase  Correlation  Representations 

The  Covariance  matrix  elements  are  directly  related  to  polarimetric  radar  measurables,  comprised  of  the 
Co/Cross-Polar  Power  Densities  Pc{p\  Px{p),  ^  (p)>  and  the  Co/Cross-Polar  Phase  Correlations 

rXp\  rXp\  rXp)  [81],  where 


[C,e  (p)] 


P,( P)  V2  R,(p)  Rc( p) 

■PrM  2P,(p)  PrM' 

rM  P  ^(p)  P/(p) 


(3.68) 


Once  the  covariance  matrix  has  been  measured  in  one  basis,  e.g.,  [C3i(//,F)]  in  {H  V}  basis,  it  can  easily 

be  determined  analytically  for  any  other  basis  by  definition  of  (3.60).  Plotting  the  mean  power  returns  and 
phase  correlations  as  function  of  the  complex  polarization  ratio  p  or  the  geometrical  polarization  ellipse 
parameters  \ff ,  X  ,  of  (3),  yields  the  familiar  ‘ polarimetric  signature  plots  In  addition,  the  expressions  for 


the  degree  of  coherence  //(/?)  and  polarization  Dp  (/?)  defined  in  (2.30)  and  (2.31),  respectively  are  given 
according  to  [34]  by 


,"(/■>)= 


M pWIp) 


Dp(p)= 


[p,(p)-pM 
(pAp)+pM 


2  \7i 


,  where  0  <  p(p) <  D  {p) <  1 


(3.69) 

and  for  coherent  (deterministic)  scatterers  p  =  Dp  =  1 ,  whereas  for  completely  depolarized  scatterers 


P  =  Dp  =  0. 


The  covariance  matrix  possesses  additional  valuable  properties  for  the  reciprocal  back-scattering  case  which 

f  .  \ 


can  be  demonstrated  by  transforming  [C3/  (//,  V )]  into  its  orthogonal  representation  for  p  = 
that 
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leading  to  the  following  inter-channel  relations 

Rt{p)\ 

and  the  symmetry  relations 
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(3.70) 


(3.71) 


(3.72) 


RTO-EN-SET-081 


1  -37 


Basics  of  SAR  Polarimetry  I 


ORGANIZATION 


Similar,  but  not  identical  relations,  could  be  established  for  the  Pauli-Coherency  Matrix  [C3P(yC>)] ,  which  are 
not  presented  here.  There  exists  another  polarimetric  covariance  matrix  representation  in  terms  of  the  so- 
called  polarimetric  inter-correlation  parameters  cr0,  p,  8,  P  ,  y ,  and  £  ,  where  according  to  [19,  Chapter 

5] 


i  p4s  pp 
p*4s  8  s  -\]y8 
p  Py  a  Py  8  Y 


(3.73) 


with  the  polarimetric  inter-correlation  parameters  cr0 ,  p ,  8 ,  /?  ,  y ,  and  £  defined  as 


(3.74) 


This  completes  the  introduction  of  the  pertinent  polarimetric  matrix  presentations,  commonly  used  in  radar 
polarimetry  and  in  polarimetric  SAR  interferometry,  where  in  addition  the  polarimetric  interference  matrices 
need  to  be  introduced  as  shown  in  [19],  after  introducing  briefly  basic  concepts  of  radar  interferometry  in 
[70]. 

3.12  Alternate  Matrix  Representations 

In  congruence  with  the  alternate  formulations  of  the  of  the  polarization  properties  of  electromagnetic  waves, 
there  also  exist  the  associated  alternate  tensorial  (matrix)  fonnulations  related  to  the  ‘four-vector 
Hamiltonian’  and  ‘ spinorial’  representations  as  pursued  by  Zhivotovsky  [109],  and  more  recently  by 
Bebbington  [32].  These  formulations  representing  most  essential  tools  for  dealing  with  the  ‘general  bi-static 
(non- symmetric)  scattering  cases’  for  both  the  coherent  (3-D  Poincare  and  3-D  Polarization  spheroid)  and 
partially  coherent  (4-D  Zhivotovsky  sphere  and  spheroid)  interactions,  are  not  further  pursued  here;  but  these 
‘more  generalized  treatments  ’  of  radar  polarimetry  deserve  our  fullest  attention. 
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